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Abstract: Animal models provide myriad benefits to both experimental and clinical research. Unfortunately, in many situations, they 
fall short of expected results or provide contradictory results. In part, this can be the result of traditional molecular biological 
approaches that are relatively inefficient in elucidating underlying molecular mechanism. To improve the efficacy of animal models, 
a technological breakthrough is required. The growing availability and application of the high-throughput methods make systematic 
comparisons between human and animal models easier to perform. In the present study, we introduce the concept of the comparative 
systems biology, which we define as “comparisons of biological systems in different states or species used to achieve an integrated 
understanding of life forms with all their characteristic complexity of interactions at multiple levels”. Furthermore, we discuss the 
applications of RNA-seq and ChIP-seq technologies to comparative systems biology between human and animal models and assess 


the potential applications for this approach in the future studies. 
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Accurately modeling the physiology and pathology 
of human systems research requires the establishment of 
a quality animal model (Alvarado & Tsonis 2006; 
Francia et al, 2011, 2006; Gotz & Lttner 2008; Hasenfuss 
1998; Lieschke & Currie, 2007). To this end, generally, 
how closely the model should mimic the human disease 
depends on the scientific question under investigation. 
Only in cases when the causal connections—structure 
function relationship or regulation of gene expression— 
are definitive, can the differences between human and 
animal models have minor effect on the analysis results 
(Hasenfuss, 1998). For example, although the zebrafish 
(Danio rerio) 1s phylogenetically distant from humans, 
its use as a complete animal model for in vivo drug 
discovery and development is growing rapidly 
(Chakraborty et al, 2009). However, if the 
pathophysiological processes are studied, especially for 
the complex diseases, then models should mimic clinical 
settings as closely as possible, otherwise the expected 
results may not be achieved or the findings of such 
studies will be of limited value. 

Accordingly, comparisons between human and 
animal models are becoming increasingly important for 
both clinical and fundamental applications (Alini et al, 
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2008; Cox et al, 2009; Fuentes et al, 2009; Huh et al, 
2010; Merchenthaler & Shughrue, 2005; Nestler & 
Hyman, 2010; Northoff, 2009). Among the available 
strategies to assess this connection, comparative systems 
biology has begun attracting special attention (Cox et al, 
2009). 

In this review, we introduce the concept of 
comparative systems biology. Next, we focus on the 
applications of next-generation sequencing methods, 
including RNA-seq and ChIP-seq, to comparative 
systems biology between human and animal models, 
before outlining some general directions of future 
developments and impacts of these types of studies. 


The rise of comparative systems biology 
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One of the greatest twentieth century achievements 
in biological research is undoubtedly the sequencing of 
different genomes. There are now complete genome 
sequences for more than 1,000 organisms (excluding 
bacteria and archaea), with more sequences being 
completed (Henkelman, 2010). Once the genome of a 
species is available, researchers are able to begin 
mapping sequences against humans and find candidate 
disease genes and build a proper disease model. However, 
the ability to fundamentally understand the genotype— 
phenotype relationship in a distinct species is often 
hindered by the inherent complexity of biological 
systems. The difference in genotype—phenotype 
relationships between human and animal models may 
originate from three sources (Figure 1): (1) functional 
divergence of genes or proteins; (2) gene deletions or 
duplications; and (3) divergent up- or down-stream 
components, out of which gene deletions or duplications 
may play the leading role (Jaillon et al, 2004). 


A Divergent function 


B- 站 ~ 国 


Gene A’ 





B Gene duplication of deletion 


Gene A’ ©- 
D 
Gene A” @-- 





Gene A 


C Divergent up- or downstrean 


Gene B 
=» 全 一 全 
x( 小 Gene A 
Gene A’ O- - } 
27 





他 


i . 
. N 
Gene C' ~ T ~ ) GeneC 


Figure 1 Mechanisms of different genotype-phenotype 

relationships between human and animal models 
The difference in genotype—phenotype relationships between human and 
animal models may originate from three sources: (A) functional divergence 
of genes or proteins; (B) gene deletions or duplications; and (C) divergent 
up- or down-stream components, out of which gene deletions or 
duplications may play the leading role. In the schematic drawing, Gene A 
and Gene A’ are orthologs while Gene A’ and Gene A” are paralogs due to 
gene duplication. 


Over the last decade, this third mechanism has 
received more attention in systems biology. The Rb 
(Retinoblastoma) gene family is a good case, because the 
members in this family are functionally conserved while 
the involved pathways are divergent between C. elegans 
and humans (van den Heuvel & Dyson, 2008). Likewise, 
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a previous study reported that over 20% of the essential 
genes for humans are non-essential for mice (Liao & 
Zhang, 2008). Consequently, traditional molecular 
biology techniques, while providing valuable insights 
into individual and/or simple genotype—phenotype 
relationship, are insufficient in deducing the complex 
phenotype-genotype relationships. Therefore, the more 
systematic methods at the systems biology level are 
necessary. 

The ultimate goal of systems biology is generating 
successful models to comprehensively describe living 
organisms. Comparative systems biology, an important 
subfield of systems biology, has no straightforward 
definition. In animal model research, the term first 
appeared in Ogawa et al’s (2008) work, reporting a 
comparative study of circadian oscillatory network 
models of Drosophila. Here, we define comparative 
systems biology as “comparisons of biological systems in 
different states or species to achieve an integrated 
understanding of life forms with all their characteristic 
complexity of interactions at multiple levels.” The 
comparison can be performed either horizontally (e.g., 
between individuals or states) or longitudinally (between 
species). The latter, which is mainly focused on human 
and animal models, is reviewed in detail here. 

Over the past decade, comparative systems biology 
has attracted widespread interest, especially for its utility 
in comparisons between human and animal models of 
complex diseases. Miller et al (2010) used a systems 
biology approach to find a number of divergent network 
modules relevant to Alzheimer disease between humans 
and mice. In a previous work, we compared humans and 
four common animal models of cardiovascular disease 
through comparative transcriptome and pathway analysis, 
revealing that a few pathways have functionally diverged 
(Zhao et al, 2012). A recent review highlighted that the 
emerging technologies in comparative systems biology 
between human and animal models offers a platform to 
systematically explore not only the molecular 
mechanism of a particular disease, thus leading to the 
identification of disease modules and pathways, but also 
the molecular relationships among distinct 
(patho)phenotypes (Barabasi et al, 2011). 

The majority of recent comparative systems biology 
studies on obtain their data through traditional high 
throughput technologies, such as microarray and ChIP- 
chip. Despite the experimental and statistical rigor as 
well as substantial insights gained through these methods, 
there has been a fundamental shift from these first- 
generation technologies (microarray and ChIP-chip) to 
next-generation sequencing (RNA-seq and ChIP-seq) 
over the last five years. We surmise that the applications 
of next-generation sequencing methods will serve a 
crucial function in the field of comparative systems 
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biology between human and animal models, offering a 
number of potential advantages. 


RNA-seq in transcriptome studies 


Previous studies demonstrated that changes in gene 
expression underlie many or even most of phenotypic 
differences between species (Marques et al, 2008; Yanai 
et al, 2004). As a result, comparative transcriptome 
analysis potentially provides information on functional 
conservation for candidate human disease genes within 
animal models. 

Initial trancriptomics studies largely relied on 
hybridization-based microarray technologies and have 
yielded valuable insights into the functional divergence 
between human and model animals (Enard et al, 2002; 
Liao & Zhang 2006). However, microarray technology 
has several limitations: over reliance upon existing 
knowledge about genome sequences; high background 
levels owing to cross-hybridization; and a limited 
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dynamic range of detection owing to both background 
and saturation of signals (Wang et al, 2009). Recent 
advances in the DNA sequencing technology have 
enabled sequencing of cDNA derived from cellular RNA 
by massively parallel sequencing strategies, a process 
termed RNA-seq (Garber et al, 2011; Mortazavi et al, 
2008). Compared with the microarray, RNA-seq has the 
advantage of allowing high-resolution characterization 
and quantification of transcriptomes with low 
background noise and the ability to distinguish different 
isoforms. 

Figure 2 shows the key procedures performed 
during RNA-seq analysis of comparative transcriptomes 
between human and animal models. The computational 
challenges in this process have been reviewed in detail 
by (Garber et al, 2011), therefore, we mainly illustrated 
the potential advantages of RNA-seq in comparative 
systems biology, including (a) comparisons between 
human and non-model animals, and (b) actual biological 
systems induced by the states of gene expression. 
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Figure 2 RNA-seq methods in comparative transcriptome analysis 
There are two strategies for sequencing animal models. If the genome was not complete or was badly annotated, the genome-independent approach should be 


used (right part). The genome-guided approach is more typical (left part). 


Though a variety of organisms have been 
genomically sequenced, the majority of these are used as 
model organisms. Since microarray relies on the genome 
information, this technique has serious limitations in 
both quantifying and comparing gene expression profiles 
from non-model animals. RNA-seq, meanwhile, can be 
applied to reconstruct the complete and high-resolution 
transcriptomes across all species. To build the 
transcriptome, several methods based on RNA-seq have 
been developed, usually falling into two main classes: 
the ‘genome-guided’ (Guttman et al, 2010; Trapnell et al, 
2010) and genome-independent classes (De novo 
assembly) (Birol et al, 2009; Schulz et al, 2012). The 
first methods rely on a reference genome to initially map 
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all the RNA-seq reads to the genome and then assemble 
overlapping reads into transcripts. Unfortunately, the 
genome-guided method is not always effective, both 
because despite a large drop in the cost of next- 
generation sequencing, the study of a complete genome 
is still costly and difficult, especially for non-model 
organisms, and because the particular model being 
studied may be sufficiently different from its reference 
genome because it comes from a different strain or line. 
Consequently, de novo assembly is particularly suitable 
for application to obtain accurate reconstructions. A 
recent study reported a large RNA-seq data set obtained 
from six organs of nine different mammals (human, 
chimpanzee, bonobo, gorilla, orangutan, macaque, 
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mouse, opossum, and platypus) and one bird (chicken), 
including both males and females (Brawand et al, 2011), 
demonstrated the utility of applying comparative systems 
biology between human and non-model animals and 
elucidated the large evolutionary gaps among these 
model organisms. 

Determining the expression states (1.e., the presence 
or absence) of genes with low abundance is a challenge 
for microarray. Consequently, the reconstruction of the 
actual biological networks (e.g., protein-protein 
interaction, transcriptional regulation network, or 
metabolic network) in either human or animal models in 
a specific condition is very difficult, not to say anything 
of the difficulty in comparing the dynamic networks 
(Farmer et al, 2012). Moreover, abnormal variations in 
alternative splicing are also implicated in disease, thus 
alternative splicing is a critical factor to consider in 
building a proper and viable animal model (Luco et al, 
2011). Unfortunately, obtain the precise alternative 
splicing map using the microarray technique 1s almost 
impossible. 

RNA-seq data is highly replicable with relatively 
little technical variation. For many purposes, RNA-seq 
may be sufficient to sequence each mRNA sample once. 
The information obtained in a single lane of RNA-seq 
data appears to be comparable to that in a single array, 
and is therefore useful in enabling the identification of 
differentially expressed genes and allowing for additional, 
further analyses, such as detection of low-expressed 
genes, novel transcripts and alternative splice variants. In 
using this method, researchers can obtain actual 
biological networks in both human and animal models, 
and garner biologically meaningful results by comparing 
between these two networks. Rowley et al (2011), for 
example, compared the actual transcriptome in platelets 
between humans and mice, providing critical information 
used in the design of mouse models of hemostasis and in 
catalyzing the discovery of new platelet functions.. 


ChIP-seq for detecting regulation changes 


Molecular interactions between proteins and DNA 
play an essential role in the regulation of gene expression 
(Cawley et al, 2004; Pokholok et al, 2006). Accordingly, 
changes in protein-DNA interactions between human 
and animal models may lead to the divergent functions of 
homologous pathways (Brown et al, 2011; Greber et al, 
2010), which is also an important aspect of comparative 
systems biology. 

Chromatin immunoprecipitation (ChIP) followed by 
genomic tiling microarray hybridization (ChIP-chip) has 
become the most widely used approach for genome-wide 
identification and characterization of in vivo protein- 
DNA interactions during the past decade (Ho et al, 2011). 
Specifically, when applied to the study of animal models 
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of human disease, CHIP-chip approaches led to many 
important discoveries in relation to transcriptional 
regulation (Chen et al, 2008), epigenetic regulation 
through histone modification (Heintzman et al, 2007), 
and evolution of protein-DNA interactions (Kim et al, 
2007). 

Like the microarray technique, CHIP-chip also has 
some limitations arising from the innate characteristics of 
microarray hybridization. Chromatin immunoprecipitation 
followed by sequencing (ChIP-seq) makes it possible to 
obtain the accurate information about the genome-wide 
profiling of DNA- protein interaction. Compared to the 
CHIP-chip, ChIP-seq has a higher resolution, fewer 
artifacts, a larger coverage and a more extensive dynamic 
range (Blow et al, 2010; Johnson et al, 2007; Mardis 
2007; Schmid & Bucher, 2007; Visel et al, 2009). 
Subsequently, we will introduce the practical 
applications of ChIP-seq in comparison between human 
and animals, including (1) identifying the regulatory 
sequences, and (2) tracing the evolution of epigenetic 
regulation. 

The human genome project, while obtaining the 
complete genomic sequences, leaves open the question 
of how to identify the regulatory sequences that control 
the spatial and temporal expression of genes unanswered 
(Birney et al, 2007; McGaughey et al, 2008). Through 
applying the ChIP-seq techniques with the enhancer- 
associated protein p300 from mouse embryonic heart 
tissue, Blow et al (2010) made an attempt to identify 
candidate heart enhancers on genomic scale, revealing 
that most of the candidate heart enhancers were less 
deeply conserved in vertebrate evolution when compared 
to the enhancers that are active in other tissues. Such 
methods could also be applied to identification of other 
transcriptional factors (TFs), and therefore are helpful in 
the reconstruction of the transcriptional regulation 
network in human and animal models. Thankfully, the 
decreasing cost of ChIP-seq has extended the 
comparative systems biology investigation to some TFs. 
For example, Schmidt et al (2010) used ChIP-seq to 
determine experimentally the genome-wide occupancy of 
two TFs, 1.e., CCAAT/enhancer-binding protein alpha 
and hepatocyte nuclear factor 4 alpha, in the livers of 
five vertebrates, revealing large interspecies differences 
in transcriptional regulation and providing insight into 
the evolution of regulatory networks. 

Epigenetic regulation is now accepted as being 
closely associated with human development, and 
subsequently many developmental disorders may be 
caused by the dysfunction of this regulation (Gottesman 
& Hanson, 2005). However, due to the deficient 
knowledge of this phenomena in other animals, build 
proper animal models for these studies is difficult. 
Nevertheless, a recent study that employed the CHIP-seq 
technique to investigate the epigenetic regulation of 
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histone H3 K4 on frogs (Xenopus tropicalis), revealed a 
hierarchy in the spatial control of zygotic gene activation 
(Akkers et al, 2009). Taken together, these advances lead 
us to speculate that the applications of CHIP-seq in 
comparative systems biology will be of great help in 
understanding embryonic diseases. 

Despite the advances that ChIP-seq offers, 
researchers should be cautious when performing ChIP- 
seq analysis because the experimental steps in ChIP-seq 
involve several potential sources of artefacts (Park, 
2009). For example, one challenge in this technique is 
that the identified enriched regions are of different types 
for different proteins (for details, refer to (Park, 2009)). 
The other potential source of artefacts comes from the 
divergence of both protein and DNA; therefore when 
using this analysis, the control experiment should be 
designed carefully. 


Perspective applications 
systems biology 


of comparative 


Comparative systems biology takes advantage of 
the systematic information from other organisms and can 
be used to great effect in studying human physiology and 
disease. Over the coming years, we expect many exciting 
developments as this field evolves in several potential 
directions. 


Dynamic networks 


Biological systems exhibit complex dynamic 
behavior, enabling cells to react to various conditions or 


cell states such as cell cycle progression (Zhu et al, 2007). 


Although static biological systems have been well 
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Abstract: In animal societies, some stressful events can lead to higher levels of physiological stress. Such stressors, like social rank, 
also predict an increased vulnerability to an array of diseases. However, the physiological relationship between social rank and stress 
varies between different species, as well as within groups of a single species. For example, dominant individuals are more socially 
stressed at times, while at other times it is the subordinate ones who experience this stress. Together, these variations make it difficult 
to assess disease vulnerability as connected to social interactions. In order to learn more about how physiological rank relationships 
vary between groups of a single species, cortisol measurements from hair samples were used to evaluate the effects of dominance 
rank on long-term stress levels in despotic and less stringent female rhesus macaque hierarchal groups. In despotic groups, cortisol 
levels were found not to be correlated with social rank, but a negative correlation was found between social rank and cortisol levels 
in less stringent hierarchies. Low ranking monkeys in less stringent groups secreted elevated levels of cortisol compared to higher 
ranking animals. These data suggest that variations in the strictness of the dominance hierarchy are determining factors in rank 
related stress physiology. The further consideration of nonhuman primate social system diversity and the linear degree of their 
hierarchies may allow for the development of valid rank-related stress models that will help increase our understanding and guide the 
development of new therapeutics for diseases related to human socioeconomic status. 
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In many animal societies, stress generated within 
long-term dominant-subordinate relationships often 
predicts an increased vulnerability to an array of diseases, 
e.g., cardiovascular disease, rheumatoid arthritis, 
respiratory, reproductive, immune related disorders and 
mental illnesses (Adler et al, 2000; Kawachi & Kennedy, 
2006; Shively & Clarkson, 1994; Siegrist & Marmot, 
2004; Wilkinson, 2001; Cohen et al, 1997; Manuck et al, 
1995; Sapolsky & Share, 1994; Ader & Cohen, 2001; 
Dhabhar & McEwen, 1999). In response to varying 
stressors, animals release hormones and glucose into the 
bloodstream to provide extra energy, which also inhibits 
digestion, growth, tissue repair and reproduction. This 
secretion is more commonly referred to as the “fight or 
flight” reaction, or the stress response, which is 
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accompanied by enhanced cognition, enhanced immune 
defenses, blunted pain perception, and sharpened sensory 
thresholds (McEwen & Lasley, 2004; Sapolsky, 2004b). 
This stress response helps animals adapt to acute stress 
and/or short-term reactions to stress. Chronic stress, 
however, such as that induced by a social hierarchy, leads 
to excessive activation of this response, which, in turn, 
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has been shown to have pathogenic ends (Sapolsky, 
2004a). 

A number of stress-related physiological studies 
have shown that animals experiencing more social stress 
due to dominance rank within a hierarchy exhibit 
hyperactivity of the hypothalamus-pituitary-adrenal (HPA) 
axis, including elevated basal levels of glucocorticoids, 
enlarged adrenal glands accompanied by increased hormone 
secretion, a repressed stress response to glucocorticoids 
when dealing with a challenge, an impaired sensitivity of 
the HPA axis to negative feedback regulation and a 
delayed recovery from stress (Sapolsky, 2005). 

The effects of dominance rank on subordinate and 
dominant individuals vary between species, and even 
within species. These variations largely depend on the 
dominance style, the means of maintaining a despotic 
dominance, the style of breeding system, rank stability, 
the availability of coping outlets for subordinates, the 
ease with which subordinates avoid dominant individuals, 
the availability of alternative strategies to overt 
competition, and even individual disposition or 
personality (Sapolsky, 2005). In some cases, dominant 
individuals are more socially stressed, such as when 
species maintain the hierarchy through frequent physical 
reassertion of dominance, in species that have a 
cooperative style of breeding system, those with transient 
periods of major rank instability, and species where 
dominants perceive neutral interactions as challenging 
while subordinates take advantage of coping strategies 
(Sapolsky, 2005). By contrast, these profiles are 
commonly observed in species with several defining 
characteristics: among subordinates in species where a 
despotic hierarchy is maintained through nonphysical 
intimidation, 1n species with consistently stable ranks, 
lower coping outlets available for subordinates, a lack 
alternative strategies to overt competition, and in species 
where dominants are adept at exerting social control and 
are highly affiliative whereas subordinates are poor at 
exploiting opportunities for social support, and likewise 
for species kept in an enclosure without sufficient space 
for subordinates to evade dominant animals (Sapolsky, 
2005). 

The two main influences of the stress response are 
the activation of the sympathetic nervous system and the 
release of glucocorticoids (such as cortisol or 
hydrocortisone in primates, and corticosterone in most 
rodent species) by the HPA axis (Sapolsky et al, 2000). 
In response to acute stressors, the sympathetic nervous 
system secrets catecholamine hormones within seconds, 
making it difficult to obtain accurate measures of the 
acute stress response. Meanwhile, glucocorticoid levels 
change over the course of 1—2 minutes, allowing for an 
extensive study of the correlation between social rank 
and glucocorticoid changes. While both glucocorticoids 
and catecholamines are essential for animals to survive 
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acute stressors, they are pathogenic when secreted in 
excess (Sapolsky et al, 2000). As such, animals of low 
social rank have been found to present a pathological 
condition consisting of elevated basal glucocorticoid 
levels, a slowed on/off “switch” of stress response, and a 
blunted sensitivity to excessive levels of glucocorticoids. 
This pathological condition has been observed in 
cynomolgus monkeys (Adams et al, 1985), talapoin 
monkeys (Keverne et al, 1982), olive baboons (Sapolsky, 
1990), squirrel monkeys (Manogue et al, 1975) and 
lemurs (Schilling & Perret, 1987). Although stress in 
subordinates has been shown to be related to the 
hypersecretion of glucocorticoids, there are exceptions 
reported among macaques (Bercovitch & Clarke, 1995; 
Chamove & Bowman, 1976; Gust et al, 1991; Gust et al, 
1993; Van Schaik et al, 1991), squirrel monkeys (Coe et 
al, 1979; Mendoza et al, 1978; Steklis et al, 1986), 
marmosets (Saltzman et al, 1994), talapoin monkeys 
(Keverne et al, 1982) and ring-tailed lemurs (Cavogelli, 
1999). 

As the primary glucocorticoids in humans, non- 
human primates and many larger mammals, cortisol has 
been used to assess stress response levels in the plasma, 
saliva, urine, and feces of many species (Cattet et al, 
2003; Constable et al, 2006; Keay et al, 2006; 
Millspaugh et al, 2002; Wallner et al, 1999; Whitten et al, 
1999; Weinegrill et al, 2004; Abbott et al, 2003). The most 
commonly used assays to detect cortisol levels in these 
samples include radioimmunoassays (RIAs), liquid 
chromatography-mass spectrometry (LC-MS/MS) and 
enzyme-linked immunosorbent assays (ELISA) (Gatti et 
al, 2009). Although plasma samples provide a 
measurement of cortisol levels at a single point in time, 
they can only be used to assess acute stress, and are 
susceptible to physiological fluctuations. Similarly, 
cortisol levels measured in the plasma peak in the early 
morning and gradually decrease to their lowest levels in 
the evening, necessitating multiple blood samples being 
taken waking until of sleep in order to accurately assess 
stress levels. Though precise, this methodology 1s difficult 
to apply to larger populations, and the compliance of 
individual participants with the requisite sampling 
schedule varies (Russell et al, 2011). Likewise, the 
cortisol measured in blood samples reflect total cortisol, 
and can thus be affected by changes in levels of cortisol- 
binding globulin (e.g., by birth control pills or 
pregnancy), which give the illusion of increased stress 
levels. Furthermore, the act of blood sampling via 
venipuncture may be a source of stress and increase 
plasma cortisol levels (Vining et al, 1983), potentially 
yielding results that do not reflect an unbiased 
observation. 

Contrary to the more invasive methods used in 
plasma samples, salivary, urinal and fecal samples can be 
obtained relatively more easily and the measured cortisol 
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reflects free cortisol. Unfortunately, the concentrations of 
cortisol still fluctuate significantly throughout the course 
of the day (Beerda et al, 1996; Cattet et al, 2003; Ekkel 
et al, 1996; Keay et al, 2006; Mormede et al, 2007; 
Owen et al, 2005; von der Ohe & Servheen, 2002). To 
overcome the issue of an animal’s diurnal rhythm, 
samples must be collected all day (Burch, 1982). More 
importantly, the cortisol levels measured in these media 
only reflect short-term stress occurring over hours to 
days, so without repeated sampling of animals, it is not 
possible to assess chronic stress levels that occur over 
weeks or months (Keay et al, 2006; Owen et al, 2005). 

Hair offers some interesting alternatives in that it 
has a fairly predictable growth rate of approximately 1 
cm/month. Likewise, hair follicles can be used to track 
cortisol levels over time, with the most proximal 1 cm 
segment to the scalp approximates the last month’s 
cortisol production, the second most proximal 1 cm 
segment approximates the cortisol production during the 
month before that, and so on (Wennig, 2000). This makes 
cortisol measured in hair a biomarker of chronic stress, 
as demonstrated previously in two studies of rhesus 
macaques in which hair samples were obtained to assess 
long-term stress levels (Davenport et al, 2006; Feng et al, 
2011). Moreover, hair samples are easily transported and 
stored in envelopes or vials at room temperature for 
years, and gathering the samples does not affect stress 
levels of the subjects (Gow et al, 2010; Russell et al, 
2011). 

In the present study, we assessed long-term stress 
levels by measuring cortisol concentrations in hair 
samples from female rhesus macaques in order to test the 
predominant hypothesis of whether stress levels may be 
related or influenced by the strictness of the dominance 
hierarchy. 


MATERIALS AND METHODS 


Subjects 

In total, 20 female rhesus macaques (Macaca 
mulatta) living in five breeding groups (n=5 in each 
group, where one monkey in the group was male) at the 
Kunming Primate Research Center, Chinese Academy of 
Sciences were used in this study. The monkeys ranged 
from 12 to 20 years of age (14.75+1.89 years), and were 
housed in colonies with access to a connected indoor 
(2.61<2.46x2.58 m)-outdoor (2.67x2.66<2.67 m) cage. 
All animals were given commercial monkey biscuits 
twice a day with tap water ad libitum, and were fed with 
fruits and vegetables once daily. All subjects had lived in 
their respective social groups at least 1 year prior to 
initial observation. All animal procedures were approved 
by the Institutional Animal Care and Use Committee of 
Kunming Institute of Zoology and were carried out in 
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accordance with the National Institute of Health Guide 
for the Care and Use of Laboratory Animals. 


Experimental design 

Animal behaviors were video recorded using a focal 
follow technique and analyzed to calculate the linearity 
of hierarchies and dominance rank for each monkey 
(Altmann, 1974). After completion of the video 
recordings, hair samples were obtained to measure the 
cortisol levels. Afterward, correlations between social 
rank and hair cortisol levels were calculated. 


Behavior sampling 

The monkeys were given 7 days to acquaint 
themselves with the observers and cameras prior to 
recording and sampling, at which time a digital camera 
fixed on a tripod was set up in front of the colony to 
record one of the monkeys in the cage. Observers kept as 
far away as possible (minimum of 5 meters) from 
enclosures to avoid disturbing the animals during 
recording. Fourteen 1-hour recordings were collected for 
each monkey throughout a 7 day period. Two recordings 
were collected per day, one taken for two hours in the 
morning (9:00-11:00) and another segment in the 
afternoon (14:00-16:00). All video recordings were 
stored on a hard disk before being reviewed and 
interpreted by three technicians. The three viewers 
analyzed each video recording simultaneously and came 
to a consensus regarding the behavior classification. 


Social rank determination 

The social ranks for monkeys in each group were 
calculated according to David’s score (DS), which is 
based on the consistent outcomes of agonistic encounters 
(Gammell et al, 2003; Zhao et al, 2011). Aggressive 
behaviors were divided into extreme aggression and mild 
ageression. Behaviors that involved physical contact (a 
bite, slap or grab)were classified as extreme aggression. 
A stare threat, open-mouth threat, chase, or displacement 
was considered a mild aggression. Submissive behaviors 
were also sub-divided into extreme and mild. Extreme 
submissive behaviors included a scream, a scream threat, 
crouch, or fleeing. Mild submissions included a lip 
smack, grimace, submissive present, or moving away 
(Shively et al, 2005). 

The averaged frequencies of agonistic behaviors 
were used to assign a hierarchal rank for each monkey 
based on their individual DS. A detailed method for 
calculating a DS has been described previously 
(Gammell et al, 2003). Standardized DS scores (DSs) 
between each group were then produced in the following 
manner: the smallest DS value (the largest minus value) 
in a group was added to its absolute value to create a 
zero value, then other DS values in the same group were 
added with the same value. The new DS values were then 
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divided by the largest new DS value. For each group, this 
process created DS; values ranging from 0 to 1, 
indicating the lowest rank (DS's = 0) to the highest rank 
(DSs = 1) in the respective groups. 

Furthermore, an H value calculated using a 
previously described equation (Singh et al, 1992; Singh 
et al, 2003) to determine the linearity of the social 
hierarchy. Calculated H values ranged from 0 to 1, 
indicating the social hierarchy on a continuum from the 
total absence of a ranking system (H=0) to a perfect 
linear order (H=1). 


Hair sampling and cortisol extraction 

Hair samples from all monkeys were collected 
between 13:30 and 15:00. At 13:30, each monkey was 
captured by an experienced technician using a net and 
removed from the colony. The hair was then taken from 
the back of the monkey’s neck using a pair of scissors, 
with particular attention made by technicians to not 
break or damage the skin. Hair was placed into a small 
pouch of aluminum foil for protection and stored as 
previously described (Davenport et al, 2006; Wennig, 
2000). 

The hair cortisol extraction was done as described 
in detail in previous studies (Davenport et al, 2006; Feng 
et al, 2011). Briefly, 500 mg of hair sample was washed 
twice for 3 min each time in 10 mL isopropanol to 
remove surface contaminants, dried at 37 °C for 8 hours, 
and then pulverized using a Retsch ball mill (Retsch 
M400) at 26 Hz for 2.5 min. After, 400 mg of the powdered 
hair was weighed and incubated in 8 mL of methanol at 
room temperature for 24 hours with a slow rotation to 
extract cortisol. Samples were then centrifuged at 8 000 
r/min for 5 minutes, and 4 mL of the supernatant was 
pipetted into a centrifuge tube and dried under a stream 
of nitrogen gas. The precipitate was reconstituted with 
0.5 mL of phosphate buffered saline solution and stored 
at —20 °C until assayed. The cortisol concentration in 
each sample was quantified with a radioimmunoassay kit 
(Cortisol RIA DSL-2000, America). The cortisol RIA 
was performed at the Radioimmuno Laboratory of the 
Second Affiliated Hospital of the Kunming Medical 
College. The cortisol extraction and RIA analysis were 
performed under a double-blind design, with each hair 
sample tested twice and the mean of the two hair cortisol 
values used to reduce measurement error. 


Data analysis 

Data analysis was conducted using SPSS (SPSS inc, 
Chicago, IL, USA). Spearman correlations were used to 
assess the effects of age on social rank and hair cortisol 
levels and evaluate the correlation between social rank 
and hair cortisol levels. In all analyses, P-values were 
determined from two-tailed tests, with the significance 
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level set at P<0.05. 


RESULTS 


Hierarchy linearity and social rank 

The linearity of the social hierarchy of five female 
rhesus macaque breeding colonies was measured using 
the H values, which were 0.86, 0.49, 0.86, 0.99 and 1.00 
among the different groups. Two of the five groups were 
found to have more stringent linear hierarchies, with H 
values above 0.90. The remaining three groups with H 
values less than 0.90 were considered to have a weak 
linear order in their hierarchy. 
Effects of age on social rank and hair cortisol levels 

In all the female monkeys (n=20), age did not affect 
the social rank as measured by DS values (7=—0.14, 
P=0.57) and DSs values (7=—0.13, P=0.58). Similarly, 
there were no significant relationships between age and 
hair cortisol levels within the female monkeys (7=—0.06, 
P=0.82). In despotic female monkeys (n=8), age did not 
affect the social rank, including DS values (7=—0.18, 
P=0.67) and DSs values (r=—0.18, P=0.67), and there 
were no significant relationships between age and hair 
cortisol levels within any female monkeys (7=0.25, 
P=0.55). Likewise, in less stringent hierarchies (n=12), 
age did not affect the social rank, including DS values 
(7=—0.13, P=0.68) and DSs values (r=—0.09, P=0.78). 
Overall, there were no significant relationships between 
age and hair cortisol levels within all the female 
monkeys included in our study (r=—0.25, P=0.44). 
Correlation between social rank and hair cortisol 
levels 

No significant correlation was discovered in the 
relationship between social rank and hair cortisol levels 
within all the female monkeys (n=20) (Figure 1A; 
r=—0.01, P=0.98; Figure 1B; ~=—0.06, P=0.79). However, 
a unique correlation between social rank and cortisol 
levels emerged after the female groups were separated 
into two groups according to their calculated H values. 
These two groups included monkeys in despotic female 
groups (n=8) and those in less stringent hierarchies 
(n=12). No significant correlation between social rank 
and cortisol levels was found in the despotic female 
groups (Figure 2A; r=0.35, P=0.40; Figure 2B; 7=0.35, 
P=0.40) but we did find a significant correlation between 
social rank and cortisol levels within the less stringent 
hierarchies (Figure 3A; 7=—0.62, P=0.03; Figure 3B; 
r=—0.61, P=0.03). 


DISCUSSION 


In the present study, we first examined female rank 
physiology (as quantified by cortisol levels in hair) in 
adult rhesus monkeys, then, after finding no relationship 
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Figure 1 Correlation between social rank and hair cortisol levels in all female monkeys (n=20) 
X-axes refer to DS value (A) and DSs value (B). Y-axis refers to the mean hair cortisol (ug/dL). 
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Figure 2 Correlation between social rank and hair cortisol levels in despotic female groups (n=8) 


X-axes refer to DS value (A) and DS's value (B). Y-axis refers to the mean hair cortisol (ug/dL). 





Mean hair cortisol (up/dL) 


-10 


—5 0 5 
DS value 


10 


Mean hair cortisol (g/dL) 


2.0°B 





0.25 


0.50 
DS value 


0.75 1.00 


Figure 3 Correlation between social rank and hair cortisol levels in less stringent female groups (n=12) 
X-axes refer to DS value (A) and DS's value (B). Y-axis refers to the mean hair cortisol (ug/dL). 


between rank and cortisol levels, the groups were 
separated into two classes according to the linear degree 
of their hierarchies: despotic groups with a linear 
hierarchy (H values>.90) and groups with less stringent 
hierarchies (H values<.90) (Li et al, 2004). After 


separation, two distinct trends emerged from the analysis. 


In less stringent groups, there was a negative correlation 
between social rank and cortisol levels, suggesting that 
low ranking monkeys were more socially stressed than 
higher ranking monkeys. Conversely, groups with 
despotic hierarchies did not display this trend. 

Overall, despotic groups were much stricter than the 
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less stringent group hierarchies, which meant that 
monkeys in the latter groups were relatively more 
assertive with dominance rank than monkeys in despotic 
groups. This finding suggests that the despotic patterns 
of female rhesus monkey culture have a differentiated 
effect on their rank physiology. Likely, less despotism of 
female rhesus monkey culture leads to an increased 
physiological stress for low ranking monkeys, which, in 
turn, increased release of cortisol by an excessive 
activation of the HPA axis. 

Age may be an important factor that affects social 
rank and cortisol of animals, but in our study, age had no 
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effect on the social rank and cortisol levels of female 
monkeys. On the whole, these findings are in line with 
hypothesis that variations in the strictness of dominance 
style are the determining factors in rank physiology 
(Sapolsky, 2005). 

A hierarchy is meaningful to the animals within it 
because it can affect rank-related physiology. The 
formation of a dominance hierarchy is not only 
predetermined by differences in the attributes of animals, 
but is also produced by the dynamics of the social 
interactions (Chase et al, 2002). Hierarchical structures 
vary among species. For example, dominance hierarchies 
in nonhuman primates can be linear (1.e., A>B>C>D) or 
circular (1.e., A>B>C>A) (Li et al, 2004). In this study, H 
values were used to calculate the linearity of hierarchies 
(Singh et al, 1992; Singh et al, 2003). Generally, H 
values ranging from 0.90 to 1.00 indicate a perfect linear 
order (Li et al, 2004). The low ranking animals in groups 
with H values less than 0.90 were found to experience 
more physiological stressors (as assessed by hair cortisol 
levels) than higher ranking animals. This phenomena 
probably occurs because low ranking females in the less 
stringent groups potentially lack predictive information 
and/or a sense of control in agonistic encounters with 
their counterparts, thereby contributing to the 
physiological stressors generated from such encounters. 
Moreover, subordinate females have relatively few 
coping outlets (e.g., being able to displace aggressive 
frustration onto other females with lower ranking) to 
alleviate the stress, and this stress may be attributable to 
the relatively unstable hierarchies (Sapolsky, 2004a). 

Although extensive research in nonhuman primates 
has shown that subordinates are more likely to exhibit 
higher levels of cortisol than dominants, including 
researches in female cynomolgus monkeys (Adams et al, 
1985), male squirrel monkeys (Manogue et al, 1975), 
male olive baboons (Sapolsky, 1990), female talapoin 
monkeys (Keverne et al, 1982), and male mouse lemurs 
(Schilling & Perret, 1987), a number of other studies 
have not found differences in cortisol levels between 
subordinates and dominants, including studies in female 
squirrel monkeys (Mendoza & Mason, 1991), male 
rhesus monkeys (Bercovitch & Clarke, 1995) and female 
talapoin monkeys (Yodyingyuad et al, 1985). In fact, 
other investigations have even reported higher cortisol 
levels in dominants than in subordinate common female 
marmosets (Saltzman et al, 1998), and in subordinate 
male and female cotton top tamarins (Snowdon et al, 
1985). 

These nonhuman primate studies mentioned above 
highlight a great inconsistency between cortisol and rank 
in a number of species. The inconsistent results are partly 
due to species variation, but also due to different 
methodologies. Previous studies utilized plasma or urine 
samples to measure cortisol levels. Blood sampling 
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requires stressful procedures such as the capture, 
restraint or venipuncture of the animals, can elevate 
circulating cortisol levels (Vining et al, 1983), somewhat 
altering the situation they were intending to measure. 
Similarly, a shortcoming of urine samples is that they are 
easily contaminated because of uncontrolled urine 
excretion. Furthermore, measured concentrations of 
cortisol in both these media are subject to circadian 
effects. More importantly, however, these cortisol levels 
only reflect short-term stress that occur over hours to 
days, and without the repeated sampling of animals, they 
cannot assess chronic stress levels that occur over weeks 
to months long periods (Keay et al, 2006; Owen et al, 
2005). Cortisol measured in hair, meanwhile, has been 
considered a biomarker of chronic stress and has proved 
useful in two studies of rhesus macaques in which hair 
samples were obtained to assess long-term cortisol and 
stress levels (Davenport et al, 2006; Feng et al, 2011). To 
our knowledge, this is the first study using cortisol 
measurements from hair sampling to evaluate the 
relationship between social rank and long-term stress in 
female rhesus monkeys, offering novel insight into long- 
term stress that arises from social. 

Though the realities of nonhuman primate social 
systems are rather simpler than those of human social 
systems, humans are similarly considered vulnerable to 
psychosocial stress stemming from their rank in the 
socioeconomic status (SES) hierarchy (Brunner, 1997; 
Sapolsky, 2005). Low SES individuals generally 
experience worse health outcomes than individuals ranked 
above them. Potential health related outcomes include an 
increased risks of cardiovascular disease, rheumatoid 
arthritis, respiratory, respiratory, reproductive, immune 
system and psychiatric diseases, as well as an overall 
increased mortality risk (Adler et al, 2000; Kawachi & 
Kennedy, 2006; Siegrist & Marmot, 2004; Wilkinson, 
2001). Similarly, long-term relationships of dominance 
and subordination stressors generate stressors associated 
with negative health sequelae in nonhuman primates, 
such as an increased cardiovascular disease risk (Manuck 
et al, 1995; Willard & Shively, 2011), depressive 
behaviors (Shively et al, 2005), HPA axis perturbations 
(Shively, 1998; Shively et al, 1997), ovarian dysfunction 
(Shively et al, 1997), reduced hippocampal volume 
(Willard & Shively, 2011), decreased dopamine-2 
receptor function (Shively, 1998), and altered 
serotonergic function (Willard & Shively, 2011). 

Taken together, researches on both human and 
nonhuman primates both support a role for physiological 
stress in a SES-related health gradient. As such, there is 
considerable interest in the development of ranking 
related stress models to aid the understanding of human 
SES-related diseases as well as guide the development of 
new therapeutics (Shively et al, 2005; Willard & Shively, 
2011). 
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Abstract: To proceed from sensation to movement, integration and transformation of information from different senses and reference 
frames are required. Several brain areas are involved in this transformation process, but previous neuroanatomical and 
neurophysiological studies have implicated the caudal area 7b as one particular component of this transformation system. In this 
study, we present the first quantitative report on the spatial coding properties of caudal area 7b. The results showed that neurons in 
this area had intermediate component characteristics in the transformation system; the area contained bimodal neurons, and neurons 
in this area encode spatial information using a hybrid reference frame. These results provide evidence that caudal area 7b may belong 


to the reference frame transformation system, thus contributing to our general understanding of the transformation system. 
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When we attempt to grasp or reach for an object, 
our motor system requires that our sensory system 
provide information about the location of the target. 
However, different sensory stimuli are encoded in 
different reference frames: visual stimuli are captured by 
the retina, auditory stimuli are received within the head, 
and tactile stimuli arise from the body. Movement must 
be directed towards the target depending on where it is 
positioned with respect to the body part. To proceed from 
sensation to movement, our brain needs to integrate 
information from the different senses and transform it 
into a non-retinocentric frame of reference, such as a 
body-part-centered frame of reference. The purpose of 
this study was to investigate whether area caudal 7b was 
a part of this transformation system. 

To construct a body-part-centered visual receptive 


field, retinal signals must be integrated with other signals. 


For example, by combining retinocentric visual signals 
with eye position signals, the head-centered visual 
receptive field of the cells can be constructed (Graziano 
et al, 1997). The arm-centered visual receptive field of 
cells requires integration of not only retinocentric visual 
and eye position signals but also head and arm position 
signals. In addition to visual signals, eye position, head 
position, and vestibular signals converge on neurons in 
area 7a (Andersen et al, 1985; Andersen et al, 1990b; 
Andersen et al, 1997). These signals are combined in a 
systematic fashion using the “gain field” mechanism, In 


Science Press 


which other signals modify the retinocentric visual signal. 
Anatomically, area 7a projects to area 7b and the 
ventral intraparietal cortex (VIP), which in turn connects 
to the ventral premotor cortex (PMv). Although other 
sensory signals modulate retinocentric visual signals, the 
visual receptive fields in area 7a are still eye-centered 
(Andersen et al, 1985). In the PMv, most neurons are 
bimodal cells, 1.e., they respond to both tactile and visual 
stimuli (Fogassi et al, 1992; Fogassi et al, 1996; 
Gentilucci et al, 1983; Graziano et al, 1994; Graziano et 
al, 1997). Most visual receptive fields encode nearby 
space in a body-part-centered reference frame (Fogassi et 
al, 1992; Fogassi et al, 1996; Gentilucci et al, 1983; 
Graziano et al, 1994; Graziano et al, 1997). Additionally, 
most neurons In the VIP are bimodal (Colby et al, 1993; 
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Duhamel et al, 1998) and employ a hybrid reference 
frame (Duhamel et al, 1997). 

On the basis of previous studies, it is known that 
area 7b receives visual and somatosensory inputs from 
cortical and subcortical brain structures (Andersen et al, 
1990a; Cavada & Goldman-Rakic, 1989a) and 1s also the 
target of cerebellar output, which suggests that it may be 
important In movement (Clower et al, 2001; Clower et al, 
2005). Damage to area 7b causes spatial deficits, such as 
neglect of contralateral peripersonal space (Matelli et al, 
1984), which suggests that this area is involved in spatial 
information processing. Area 7b contains cells with 
visual and tactile properties that resemble those of PMv 
cells (Dong et al, 1994; Leinonen et al, 1979). These 
results suggest that area 7b may be one component of the 
receptive field transformation system. Do neurons in this 
area have the characteristics of the transformation system? 
Does the area employ a retinocentric reference frame 
similar to area 7a, a hybrid system similar to VIP, or 
some other reference frame? To date, no quantitative 
studies on these issues have been conducted. 

In the present study, we first qualitatively tested the 
sensory properties of neurons in the caudal portion of 
area 7b. If a neuron showed a clear visual response, the 
reference frame of the neuron was tested quantitatively. 
We plotted the visual receptive fields of neurons when 
the eyes were held still. The receptive fields were then 
replotted when the eyes were in a new position. In this 
way, we could determine whether the visual receptive 
fields were centered on the eyes. 


MATERIALS AND METHODS 


Ethics statement 

All animal care and surgical procedures were 
performed in accordance with the National Institutes of 
Health Guide for the Care and Use of Laboratory 
Animals and approved by the Institutional Animal Care 
and Use Committee at the Kunming Institute of Zoology. 

Two adult male Macaca fascicularis (4.2-6.0 kg) 
were used in the present research. They were single- 
housed under standard conditions (a 12-h light/dark cycle 
with lights on from 07:00 to 19:00, humidity at 60%, and 
temperature at 21 + 2 C in the animal colony. Their 
body weights were closely monitored during the 
experiments. 


Surgical procedures 

All surgeries were performed under strict aseptic 
conditions and, unless otherwise indicated, under ketamine 
hydrochloride (10 mg/kg) and acepromazine (0.4 mg/kg) 
anesthesia. Under deep anesthesia, one stainless steel 
recording chamber for chronic single-unit recording, 2.5 
cm in diameter, was embedded in dental acrylic over the 
right parietal lobe of each monkey. A steel bolt for 
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holding the monkeys’ head stationary during experiments 

was also embedded in the acrylic. During the same surgery 

session, one scleral search coil (Fuchs & Robinson, 1966) 
was implanted to record eye position signals. 

Post-surgical recovery lasted approximately one 
week. Both monkeys were given two additional weeks 
before behavioral training started to allow the skull to 
grow tightly around the skull screws. 

After the behavioral training was completed, the 
recording chamber was opened under aseptic conditions 
and a ~4 mm hole was drilled through the dental acrylic 
and bone to expose the dura mater. All data were 
collected from this region. 


Behavioral training 

Before experimental task training began, each 
monkey’s ad libitum daily water intake was measured. 
Based on this measurement, the monkeys were placed on 
a restricted water schedule whereby each monkey 
received most of its water as a reward for performing the 
task. If the monkey received less water than its normal 
daily amount during the experimental session, a post- 
session supplement was given. Each monkey also 
received ad libitum water for two consecutive days every 
weekend during which training and recording sessions 
were not conducted. Each monkey’s body weight was 
closely monitored during the time it was on the restricted 
water schedule. 

At the beginning of the task training, the monkey 
sat in the chair with its head and left arm fixed, and its 
eye movements were constantly monitored using a 
standard eye coil technique (Fuchs & Robinson, 1966; 
Graziano et al, 1997). As shown in Figure 1, three LEDs 
were placed 57 cm in front of the monkey at eye level, 
15 visual degrees apart. First, the monkeys were trained 
to perform the fixation task without a peripheral visual 
stimulus. The central LED was set to blink at a frequency 
of 6 Hz. If the monkey fixated on the central LED within 
a specified round spatial window (fixation window), the 
LED stopped blinking and remained lit. If the monkey 
maintained fixation for the remainder of the trial, the 
LED was turned off and the monkey received a drop of 
juice as a reward. A 3—10-second intertrial interval (ITI) 
commenced at the same time that the reward was given. 
If the monkey broke fixation at any time during the trial, 
the LED was turned off, the ITI was restarted, and the 
monkey was not rewarded. As the training progressed, 
the size of the spatial fixation window was decreased and 
the fixation period was increased. At the final stage, the 
monkeys were able to fixate within a round window of 3 
visual degrees in diameter for a 4-second fixation period. 
After they learned to fixate on the central LED, the 
monkeys were trained to fixate on the two lateral LEDs 
using the same procedure. 


Volume 34 Issue E2 


E52 JIANG, et al. 


Having learned the fixation task, the monkeys were 
next trained to perform the same task when a peripheral 
visual stimulus was presented (Figure 1 and Figure 2). 


Recording procedure 

The recording sessions started after the monkey was 
able to perform the task for more than 700 trials in one 
day. During the daily recording sessions, a hydraulic 
microdriver was mounted to the top of the recording 
chamber. A steel guide cannula (22-gauge syringe needle) 
was slowly lowered through the small hole in the 
chamber until its tip was close to the surface of the dura 
mater. Then, a varnish-coated tungsten microelectrode 
(Frederick Haer, impedance 0.8—1.0 MQ) was advanced 
slowly from the guide cannula through the dura mater 
and into the cortex. After the electrode penetrated the 
surface of the cortex, as indicated by the appearance of 
background neuronal signals, the electrode was 
withdrawn 1.5 mm, and the brain was allowed to recover 
from the indentation caused by the advance of the 
electrode for about one hour. The electrode was then 
slowly moved up and down to isolate single-unit activity. 
As the experiment progressed, the cannula and a new 
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electrode were inserted into a slightly different location 
in the small hole. 


Tests 

When a single neuron was isolated, its activity was 
constantly monitored with an oscilloscope and a 
loudspeaker. Each neuron was first tested to determine if 
it had tactile, visual, and/or auditory responses. The 
tactile response was tested by manual palpation, gentle 
pressure, and stroking with a cotton swab. Tactile 
responses on the face were studied in the dark or when 
the eyes were covered. Visual responses were studied 
with hand-held objects, such as cotton swabs, a bunch of 
plastic grapes, monkey dolls, or a metal rod. The visual 
stimulus was moved into the visual receptive field from 
above, below, or the side. The receptive field was 
measured with the most effective stimulus when the 
response became stable and clear. By presenting stimuli 
at different distances from the monkey, we also tested the 
depth extension of visual RFs. 

If a cell showed clear visual responses, it was then 
tested quantitatively with one of the following two tests: 
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Figure 1 Two experimental tests were used in the present study. A: robot arm test. In each trial, the monkey fixated on one of the 
three LEDs (FIX A, B, C), and the visual stimulus was advanced along one of the five trajectories indicated by the arrows. 
B: plastic rod test. In each trial, the animal fixated on one of the three LEDs (A, B, C). During the fixation, one of the five 
plastic rods indicated by the filled circles (stimulus positions 3 to 7) descended into view from behind a panel. The circles 


on the monkey’s head show the recording locations. 
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Figure 2 Time course of the two experimental tests used in the 
present study. The change in the vertical position of 
the stimulus time course indicates stimulus movement. 
The dotted line indicates that the duration was 
random. FIX = fixation light; STIM = stimulus. 


Robot arm test 

Only monkey 1 was tested with this type of stimulus, 
as shown in Figure 1A. At 0.5 seconds after the onset of 
fixation, a visual stimulus (a piece of blue Styrofoam, 2 
cm in diameter) was advanced towards the animal by a 
mechanized track at a speed of 8 cm/sec along one of 
five equally spaced parallel trajectories at the same level, 
15 visual degrees apart. The trajectories were 12 cm long 
and their locations were adjusted to match the visual 
receptive field being studied. The stimulus stayed at the 
end of the trajectory for 0.2 seconds before the ITI 
started. The stimulus was moved to the starting position 
of the next trajectory during the ITI. The time course of 
this test 1s demonstrated in Figure 2. The combination of 
three eye fixation positions and five stimulus trajectories 
resulted in fifteen conditions. These fifteen conditions 
were tested in an interleaved fashion, and each condition 
typically had ten trials. 

However, there were some problems with the 
stimulus. Firstly, in the robot arm test, the Styrofoam 
stimulus was always in view and some area 7b neurons 
quickly habituated to it. Secondly, the five trajectories 
used in the robot arm test were parallel to each other and 
moved straight towards the monkey. When the stimulus 
was advanced towards the monkey along a trajectory 
located peripheral to the visual receptive field, it 
produced a much larger visual angle compared with other 
trajectories in the middle of the visual receptive field 
before reaching its end point. When a larger visual angle 
was covered, a larger area on the retina was stimulated, 
which made it more difficult to directly compare results 


Kunming Institute of Zoology (CAS), China Zoological Society 


Visuospatial properties of caudal area 7b in Macaca fascicularis E53 


from different trajectories. To address these problems, we 
designed a plastic rod test in which the stimulus suddenly 
came into view and the visual angle of visual stimulus 
was kept constant. 


Plastic rod test 

This test was used on both monkeys. An array of 
five small plastic rods, arranged 15 degrees apart along a 
semicircle and controlled by five different stepping 
motors, served as visual stimuli. The center of the 
semicircle was located at the midline of the monkey’s 
face, and its radius was 5 cm. There were a total of ten 
possible locations on the semicircle (labeled 1 to 10) for 
the five rods so that the locations of the stimuli could be 
adjusted according to the location of the visual receptive 
field being tested. 

Each trial began when the animal fixated on one of 
the three LEDs. After an initial fixation period of 1.25— 
2.0 seconds for monkey 1 and 2.0-2.5 seconds for 
monkey 2, one of the five rods descended at a speed of 
10 cm/sec for 6.0 cm, emerging from behind a head-level 
panel into the monkey’s view. The rod remained 
stationary for 1.0 second in front of the animal. The trial 
was then ended and the animal received 0.2 ml of juice 
as a reward. A variable ITI began (4.0-7.0 seconds for 
monkey 1 and 2.0-3.0 seconds for monkey 2), and the 
rod was withdrawn to behind the panel. Figure 2 shows 
the time course for the test. The combination of three eye 
fixation positions and five stimuli (rods) yielded fifteen 
conditions. These fifteen conditions were run in an 
interleaved fashion, and each condition typically 
contained ten trials. 


Statistical procedures 

Two analysis windows were chosen for each cell. 
The first, termed the "baseline window" and used to 
calculate the pre-stimulus baseline activity, consisted of 
the time period from 0.6 to 0.0 s before the onset of the 
stimulus. The other window, termed the "response 
window", was used to examine the neuron’s response to 
visual stimuli and consisted of the time period from the 
onset of the neuron’s main visual response (1.e., three 
standard deviations above the baseline for three 
consecutive bins, with a 50-ms bin length) to 1.5 s after 
stimulus onset because all recorded neurons continued 
their responses to the end of the trial. We chose not to use 
a fixed response window as in previous studies because 
the size and location of the receptive fields differed 
among the neurons, and as a result, the exact time when 
the stimulus moved into the receptive field varied from 
neuron to neuron. Because cells did not respond to the 
stimulus at the same time across different locations, we 
chose the response window from the earliest onset 
(which was usually also the best response) location and 
applied this window to all fifteen experimental 
conditions. The response window was shortened and 
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lengthened for several tested neurons; very similar 
results were obtained, which demonstrated the robustness 
of this method of analysis. 


Qualitative analysis 

All neurons were classified qualitatively. Firstly, 
neurons without a clear peak (as judged by the existence 
of a response location in the curve that was significantly 
different from the response for at least one location on 
both sides of the peak, according to a t-test) in three 
tuning curves were not included in the data set for the 
spatial shift analysis. Secondly, before calculating the 
shift index, we compared the shapes of the tuning curves 
by moving two of the tuning curves until the peaks of 
three tuning curves overlapped, because shift index 
values can be affected by shapes of the tuning curves. We 
classified neurons into the following three types based on 
this comparison: (1) Neurons without obvious shape 
changes (as judged by the overlap of the standard 
deviations at each stimulus location) and receptive fields 
that did not shift with the eyes were classified as non- 
retinocentric neurons. (2) Neurons without obvious shape 
changes and receptive fields that shifted with the eyes 
were defined as retinocentric neurons. (3) Neurons with a 
partial spatial shift, which included two groups of 
neurons: the first group had no obvious shape changes 
but had partially shifted receptive fields; the second 
group had obvious shape changes and was also classified 
as partial spatial shift neurons. 


Quantitative analysis 

Based on the qualitative analysis, neurons with clear 
peaked responses were analyzed quantitatively. A shift 
index (SI) was calculated to measure how much the 
visual receptive field shifted when eye position changed 
as the monkey fixated on different LEDs (Graziano, 
1999; Schlack et al, 2005). 

First, a center of mass (CM) was calculated for each 
of the tuning curves corresponding to the right (CMR) 
and left (CML) fixation LEDs. When the animal fixated 
on the right LED, the responses (number of spikes per 
second) from each of the five stimulus locations were 
numbered from the monkey’s left to right as R1 to RS. 
The formula used to calculate the CMR of the right 
fixation LED tuning curve is CMR=(R1+2xR2+ 
3x R3+4xR4+Sx R5)/(RI+ R2+R3 4+ R4+R5). 

The CML of the left LED could be calculated In the 
same way. From the CMs, the shift index (SI) was 
calculated as ST = CMR— CML. 

A positive SI indicated that the visual receptive field 
shifted in the same direction as the change in eye 
position in the orbit. A negative SI indicated that the 
visual receptive field was displaced in the opposite 
direction from the change in eye position. An SI value 
close to zero indicated a stationary visual receptive field. 
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Histology 

After the series of experiments was completed, 
monkey | was injected with an overdose of pentobarbital 
sodium (100 mg/kg) and perfused with saline and 4% 
paraformaldehyde. The head was removed from the body 
and placed in a stereotactic apparatus. The skull and dura 
mater were removed to expose the brain. The 
intraparietal, superior temporal and lateral sulci were 
mapped stereotactically. Figure 3 shows the entry 
location of the electrode in relation to the sulci for 
monkey 1. On this basis, we confirmed that the recording 
sites were located in the caudal portion of area 7b. 


MONKEY 


SPL 





è Location of electrode penetration 


Figure 3 Locations of area 7b and electrode penetration. Black 
dot: the recording area was ~4 mm in diameter and 
~2 mm in depth. IPL, the inferior parietal lobe; SPL, 
the superior parietal lobe; ips, the intraparietal sulcus. 


Monkey 2 is currently still in use in other 
experiments. Therefore, anatomical data were not 
available for this monkey. 


RESULTS 


Response categories 

In total, 291 neurons in caudal area 7b were 
examined in the right hemispheres of two awake 
monkeys. In these cells, neuronal responses were 
classified as somatosensory, visual, or bimodal 
(somatosensory + visual). We also tested 82 cells in this 
area using auditory stimuli and found that only one (1%) 
responded. Table 1 shows the proportions of the different 
response types. 


Table 1 Categories of neurons based on their response type 
Somatosensory only 15% (44/291) 
Visual only 14% (42/291) 
42% (122/291) 
29% (83/291) 


1% (1/82) 


Somatosensory + visual 
Unresponsive 
Auditory 


# Values in parentheses are the numbers of neurons. 
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Somatosensory cells 
To classify the tactile receptive fields, we arbitrarily 
divided the monkey body into four parts: head, arms 


(including the hands), trunk, and legs (including the feet). 


The tactile receptive fields were classified into three 
types according to the number of body parts 
encompassed by the receptive fields. The tactile 
receptive fields were also classified into three types: 
ipsilateral, contralateral, and bilateral according to the 
recording hemisphere. The distribution of the receptive 
fields of the somatosensory cells is shown In Table 2. 


Table 2 The distribution of tactile receptive fields 


Bimodal cells Somatosensory cells 


One body part 39% (46/119) 73% (32/44) 
Two body parts 44% (52/119) 23% (10/44) 
Three or four body parts 17% (21/119) 4% (2/44) 
Bilateral side 70% (83/119) 39% (17/44) 
Ipsilateral side 9% (11/119) 16% (7/44) 
Contralateral side 21% (25/119) 45% (20/44) 


# Values in parentheses are the numbers of neurons. 


The size of the receptive fields of the 44 
somatosensory cells varied from covering only a portion 
of one of the eyebrows to encompassing the entire body. 


Visual cells 

In total 35 of the 42 purely visual cells were 
recorded for a time to test the location of the visual 
receptive fields, whereas the remaining seven cells were 
only studied long enough to determine that they had 
visual responses. The visual receptive fields were divided 
in two ways. According to recording hemisphere, the 
visual receptive fields could be divided into three 
categories: ipsilateral, contralateral, and bilateral to the 
recording hemisphere. By presenting stimuli at different 
distances from the monkey, neurons were classified into 
two types: distance > 1 m and distance < Im. Only five 
visual neurons were tested with stimuli at different 
distances. Table 3 shows the results. 


Table 3 The distribution of visual receptive fields 


Bimodal cells Visual cells 


Bilateral side 70% (83/119) 51% (18/35) 
Ipsilateral side 8% (10/119) 11% (4/35) 
Contralateral side 22% (26/119) 37% (13/35) 
Depth > 1 m 64% ( 29/45) 40% (2/5) 
Depth < 1 m 36% (16/45) 60% (3/5) 


# Values in parentheses are the numbers of neurons. 


Bimodal cell 

We assessed the locations of the tactile and visual 
receptive fields for 119 of the 122 bimodal cells. The 
remaining three cells were only recorded long enough to 
demonstrate the existence of visual and tactile responses. 
The tactile receptive and visual receptive fields of the 
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bimodal cells were classified into different types in the 
same way as the somatosensory and visual cells, and the 
results are shown in Tables 2 and 3, respectively. 

One interesting finding is that the visual receptive 
field in bimodal neurons always surrounds the tactile 
receptive field. Overall, 76% (90/119) of the bimodal 
cells had tactile receptive fields that spatially match their 
visual receptive fields. 

We also tested the visual response to perceptual 
depth in 45 out of the 119 bimodal cells. More than half 
(64%) of these cells responded clearly to the stimuli 
presented more than 1 m away from the monkey’s face, 
and the remaining 36% of the visual receptive fields 
were confined to a depth of 1 m or less, with a mean 
inner-limit distance of 0.4 m. 

In summary, the qualitative tests showed that there 
were three major types of responsive neurons in caudal 
area 7b: visual neurons (14%), somatosensory neurons 
(15%), and bimodal neurons (42%). The receptive fields 
of these three categories of neurons were usually large: 
52% (85/163) of all tactile receptive fields (including the 
somatosensory cells and bimodal cells) encompassed two 
or more parts of the body surface, 61% (100/163) of all 
tactile receptive fields were bilateral, and 66% (101/154) 
of the visual receptive fields (including the visual cells 
and bimodal cells) were bilateral. Most visual receptive 
fields extended more than 1 m away from the body part. 
Most of the bimodal cells had spatially matched tactile 
and visual receptive fields. 


Reference frames for encoding visual stimuli 

If a neuron showed a clear visual response, it was 
further studied to determine whether the receptive field 
moved with the eye. As described in the methods section, 
two types of stimuli were used and three different types 
of neurons were found. 

For some neurons, such as Neuron JC89 (Figure 4 
and Figure 5) tested with the plastic rod stimulus, the 
visual receptive fields moved with the eye. As seen in the 
Figures, when the monkey shifted its fixation point from 
the right to the middle LED and then to the left LED, the 
position at which neuron JC89 had the best response 
shifted from position 3 to position 2 and then to position 
1. These results suggest that the visual receptive field 
moves with the eyes. 

Some neurons had visual receptive fields that did 
not shift with eye movement. Neuron HB16 (Figure 6 
and Figure 7), which was tested with the robot arm 
stimulus, belonged to this type of neuron. Figure 7 shows 
that regardless of the LED on which the monkey was 
fixating, the position with the best response did not 
change for this neuron; it always showed the strongest 
response to a stimulus presented at position 4, and the 
three tuning curves were similar and overlapping. These 
results demonstrate that the visual receptive fields do not 
move with the eyes. 
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Figure 4 Histograms of activity from neuron JC89 in area 7b. The upper row of the histograms shows the visual responses to the 
stimuli presented at the five stimulus positions when the monkey was fixating on the left LED (marked as L). The middle 
row of the histograms shows the visual responses when the monkey was fixating on the center LED (marked as M). The 
lower row of the histograms shows the visual responses when the monkey was fixating on the right LED (marked as R). 
For each histogram, the x-axis shows the time in ms, and the y-axis shows the neuronal responses in spikes/second. The 
vertical line inside each histogram indicates the time at which the visual stimulus began to move into view. The horizontal 
line inside histogram L1 indicates the temporal window used to calculate the neuronal response to the stimulus. 
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Figure 5 Three tuning curves for the activity of neuron JC89. 
The x-axis shows the stimulus position, and the y-axis 
shows the neuronal responses in spikes/second. The 
dashed line shows the average pre-stimulus activity. 
During right fixation (circles), the cell responded 
most strongly to the stimulus presented at position 3. 
During middle fixation (filled diamonds), the cell 
responded most strongly to the stimulus presented at 
position 2. During left fixation (stars), the cell 
responded most strongly to the stimulus presented at 
position 1. This result indicates that the receptive 
field of this cell shifted with changes in eye position. 
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Neuron HB86 (Figure 8, Figure 9), which was 
tested with the plastic rod stimulus, is an example of a 
neuron with properties intermediate to those of a visual 
receptive field, which did not move with the eyes, and 
those of a visual receptive field, which did. That is, these 
neurons had visual receptive fields that moved slightly. 
As seen from the three tuning curves of neuron HB86 in 
Figure 9, when the monkey’s eyes shifted from the right 
LED to the middle LED, the two tuning curves 
overlapped, 1.e., the visual receptive field did not shift 
with the eyes. In contrast, when the eye moved from the 
middle LED to the left LED, the visual receptive field 
partially shifted to the left as indicated by increased 
responses to stimuli at positions 1 and 2. 

To quantify the results, a shift index (SI) was 
calculated for each neuron to measure how much the 
visual receptive field shifted when the eye position 
changed because the monkey was fixating on different 
LEDs. An example SI calculation for neuron JC89 is as 
follows: 


CMR = (R1+2x R24+3x R34+4xR44+5x R5)/ 
(R1+R2+ R3+ R44 R5) = 3.0 


and 
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Figure 6 Histograms of the activity of neuron HB16. Axes are the same as in Figure 4. The horizontal line inside histogram L4 


indicates the temporal window used to calculate the neuronal responses to the stimulus. 


15 HB16 


Spikes/sec 





Stimulus position 


Figure 7 Three tuning curves for the activity of neuron HB16. 
Axes are the same as in Figure 5. During right LED 
fixation (circles), middle LED fixation (filled 
diamonds), and left LED fixation (stars) the cell 
responded most strongly to the stimulus presented at 
position 4. This result indicates that the receptive 
field of neuron HB 16 did not shift with changes in 
eye position. 


CML = (L1+2xL2+3xL3+4xL4+5xL5)/ 
(11+ 224+ 1234+ 14+ 25) =2.2 
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Thus, the shift index of JC89 is S/(JC89) = CMR — 
CML =3.0-2.2=0.8. 

Using the same formula, the shift indices for the 
other two example neurons are as follows: neuron HB16 
= 0.1 and neuron HB86 = 0.3. 

Overall, 34 neurons showing clear responses to the 
automatically presented stimuli demonstrated peak 
responses. We calculated the shift index for each of the 
34 neurons. Figure 10 shows the distribution of the shift 
indices was from —0.4 to 0.8. The distribution of the shift 
index was significantly greater than zero (x = 9.23, P < 
0.01). On average, the visual receptive fields moved 
when the eyes moved. Thus, some neurons had visual 
receptive fields that moved with eye movement, some 
neurons had partially moving visual receptive fields, and 
some had receptive fields that did not move with the eyes. 


Effect of eye position on the neuronal firing rate 

In addition to the location of the visual receptive 
field, the visual response magnitude of caudal area 7b 
neurons was sometimes significantly influenced by the 
position of the eyes in the orbit, which is termed the 
modulation effect. For example, regardless of the LED 
on which the monkey fixated, HB169 responded with 
almost equal strength at each of the five positions, 
although its response was slightly stronger when the 
stimulus was at position 5 (Figure 11). However, the 
response magnitude of the cell decreased monotonically 
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Figure 8 Histograms of the activity of neuron HB86. Axes are the same as in Figure 4. The horizontal line inside histogram L2 
indicates the temporal window used to calculate the neuronal response to the stimulus. 
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Figure 9 Three tuning curves for the activity of neuron HB86. 
Axes are the same as in Figure 5. During right LED 
fixation (circles) and middle LED fixation (filled 
diamonds), the two corresponding tuning curves 
overlapped. The cell showed the weakest response to 
the stimulus presented at position 1 and the strongest 
response to the stimulus at position 5. In contrast, 
when the eyes moved from the middle LED to the 
left LED (stars), the visual receptive field partially 
shifted to the left as indicated by the increased 
responses to stimuli at positions 1 and 2. Overall, the 
visual receptive field of neuron HB 86 partially 
shifted with changes in eye position. 
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Figure 10 Distribution of shift indices for 34 neurons. The 
distribution of shift indices was significantly 
above zero (P< 0.01) 


as the eye moved from the right LED to the middle LED 
and from the middle LED to the left LED. A 5 x 3 
within-subjects two-way ANOVA (5 stimulus positions x 
3 eye positions) showed that the modulation of the cell 
response was significant (P < 0.001). Neuron JC 89 
showed the same type of response modulation (P < 0.05). 

In contrast, neuron HBI6 responses were not 
modulated by eye position (P > 0.5). This cell responded 
best to the stimulus presented at position 4. The 
magnitude of the visual response did not change 
significantly when the eyes were at different positions in 
the orbit. 
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Figure 11 Three tuning curves for the activity of neuron 
HB169. Axes are the same as in Figure 5. Regardless of the 
LED on which the monkey fixated, neuron HB169 responded 
to the stimulus with almost equal strength when presented at 
each of the five positions, although the response was slightly 
stronger at position 5. However, the response magnitude 
decreased monotonically as the eyes moved from the right LED 
(circles) to the middle LED (filled diamonds) and from the 
middle LED to the left LED (stars). This result indicates that 
the visual response magnitude of neuron HB169 was 
significantly modulated by changes in eye position (P < 0.001). 


The ANOVA for the 34 neurons showed that the 
visual response magnitudes of 12 cells (35%) were 
significantly modulated by eye position, and the response 
magnitudes of the remaining 22 cells (65%) were not 
significantly modulated by eye position. 


DISCUSSION 


Visual and tactile responses in caudal area 7b 

In this study, we qualitatively tested the sensory 
properties of neurons in caudal area 7b. Three major 
types of responsive neurons were found: visual neurons 
(14%), somatosensory neurons (15%), and bimodal 
neurons (42%). The bimodal, visual-somatosensory cells 
usually responded to visual stimuli positioned close to 
the tactile receptive field. Some cells responded only to 
stimuli within centimeters of the body whereas others 
responded to stimuli presented more than 1 m away. 

Other studies regarding the general properties of 
neurons in area 7b have also demonstrated that it 
contains visual cells, somatosensory cells, and bimodal 
cells (Dong et al, 1994; Leinonen et al, 1979). The 
proportions of cells across these three categories varied 
from study to study, with that the percentage of bimodal 
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neurons varying from 6% to 33%. The present study 
showed 42% bimodal neurons. 

One possible reason for these quantitative 
differences is that the stimuli and the way in which 
stimuli are used may vary across labs. For example, we 
found that the visual response to a stimulus delivered 
from behind the animal and then moved into its view was 
usually better than the visual response observed when the 
same stimulus was delivered from in front of the animal. 
It is possible that the stimulus coming from behind was 
unexpected, and therefore, lt drew more attention. Other 
labs using stimuli presented from the front may therefore 
report a lower percentage of visual responses. Another 
possible reason for the quantitative differences is that the 
recording sites (1.e., the sample population) may differ 
from the caudal part of area 7b used in the present study. 
Furthermore, the monkeys used in our study were awake, 
whereas other studies were conducted on anesthetized 
monkeys. As described previously, the wakefulness of 
animals can affect the response properties of area 7b 
neurons (Robinson & Burton, 1980). When the monkey 
is asleep, the receptive fields of some somatosensory 
neurons become smaller, the borders become more 
difficult to define, and the somatic activations become 
less clear. 


Transformation of reference frames 

It has been suggested that body-part-centered 
receptive fields provide a general solution to the central 
problem of sensory-motor integration (Graziano et al, 
1994) and that this type of reference frame can encode 
the distance and direction from a body part to a nearby 
visual stimulus (Graziano et al, 1997). However, to 
proceed from a retinocentric reference frame to a body- 
part-centered reference frame requires a complicated 
transformation system that involves several brain areas. 

Only one study has reported that, In some cases, the 
visual responses of neurons in area 7b do not depend on 
the direction of the gaze (Leinonen et al, 1979). Our 
present research is the first quantitative study to explore 
whether the visual receptive fields of neurons in caudal 
area 7b move with the eyes. We found that neurons in 
caudal area 7b encoded space in a hybrid reference frame, 
and some neurons encoded space in a retinocentric 
reference frame. We did not further differentiate the 
reference frames employed by the non-retinocentric 
neurons, e.g., head-centered, body-centered, or world- 
centered, but showed that the caudal part of area 7b was 
one component of the reference frame transformation 
system. 

Recent studies on the inferior parietal lobe (IPL) 
challenge the widely used bipartite subdivision of the 
IPL convexity into caudal and rostral areas (7a and 7b, 
respectively). These studies provide strong support for a 
subdivision of the IPL convexity into four distinct areas 
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(Opt, PG, PFG, and PF) and indicate that these four areas 
are not only cytoarchitectonically distinct but also 
distinct in connectivity (Gregoriou et al, 2006; Rozzi et 
al, 2006). In our study, the recording site was located in 
the caudal portion of area 7b, corresponding to area PFG. 
The present results provide evidence that caudal area 7b 
may be one component of the reference frame 
transformation, which is consistent with the proposed 
function based on the connectivity pattern of area PFG 
(Borra et al, 2008; Fogassl & Luppino, 2005; Rozzi et al, 
2006). 


Eye-position-independent visual receptive fields in 
other brain areas 

The anterior intraparietal (AIP) area is located 
rostrally in the lateral bank of the intraparietal sulcus and 
is known to play a crucial role in visuomotor 
transformations for grasping. Anatomically, the AIP is 
connected to the IPL and the rostral ventral premotor 
area FS (Borra et al, 2008; Luppino et al, 1999). 
Reversible lesion experiments have shown that AIP 
inactivation produces grasping impairment (Gallese et al, 
1994). Visual, motor, and visuomotor neurons have been 
recorded during the execution of grasping movements 
(Murata et al, 1996; Murata et al, 2000; Sakata et al, 
1995; Taira et al, 1990). 

Another parietal area suggested to have eye- 
position-independent visual receptive field is the parieto- 
occipital visual area (V6A). Studies have shown that 
although most of the visual receptive fields are eye- 
centered, there is a small portion of neurons in this area 
with visual receptive fields that remain static when the 
eyes move and another small group of visual receptive 
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fields that partially move with the eyes (Fattori et al, 
1992; Galletti et al, 1993; Galletti et al, 1995). 

Area 7b, VIP, and PMv are monosynaptically 
interconnected and all project to the putamen (Cavada & 
Goldman-Rakic, 1989a; Cavada & Goldman-Rakic, 
1989b; Cavada & Goldman-Rakic, 1991; Matelli et al, 
1986). Studies of the putamen of anesthetized and awake 
macaque monkeys have shown that about 30% of the 
cells with tactile responses also respond to visual stimuli. 
For these bimodal neurons, the visual receptive fields are 
anchored to the tactile receptive fields, moving as the 
tactile receptive fields move (Graziano & Gross, 1993). 

Comparing these brain areas quantitatively by 
calculating the fraction of each cell type in each 
individual area, might provide more information about 
how the system works, and more experiments should be 
conducted in which all of these areas are tested using 
identical procedures. 

In summary, we tested the response properties and 
the spatial encoding coordinates of neurons in caudal 
area 7b. Three major kinds of responsive neurons were 
found in this area: visual neurons, somatosensory 
neurons, and bimodal neurons. The receptive fields were 
usually large, and visual receptive fields of the bimodal 
neurons were always located near tactile receptive fields. 
Neurons in this area encoded space in a hybrid reference 
frame, which may be a characteristic of the reference 
frame transformation system. In addition, the responses 
of some neurons were modulated by the eye position. 
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Abstract: The Chinese tree shrew (Tupaia belangeri chinensis) is a small experimental animal with a close affinity to primates. This 
species has long been proposed to be an alternative experimental animal to primates in biomedical research. Despite decades of study, 
there is no pure breed for this animal, and the overall genetic diversity of wild tree shrews remains largely unknown. In order to 
obtain a set of genetic markers for evaluating the genetic diversity of tree shrew wild populations and tracing the lineages in 
inbreeding populations, we developed 12 polymorphic microsatellite markers from the genomic DNA of the tree shrew. An analysis 
of a wild population of 117 individuals collected from the suburb of Kunming, China, showed that these loci exhibited a highly 
expected heterozygosity (0.616). These 12 microsatellites were sufficient for individual identification and parentage analysis. The 


microsatellite markers developed in this study will be of use in evaluating genetic diversity and lineage tracing for the tree shrew. 


Keywords: Chinese tree shrew; Microsatellite; Heterozygosity; Individual identification; Parentage testing 


The Chinese tree shrew (Tupaia belangeri chinensis) 
is a squirrel-like animal. It belongs to the Family 
Tupalldae of Scandentia. Tree shrews have the highest 
brain-to-body mass ratio of known mammals (Peng et al, 
1991). Because tree shrews share some characteristics 
with primates and insectivores, the exact phylogenetic 
position of the tree shrew has been debated (Arnason et 
al, 2002; Janecka et al, 2007; Peng et al, 1991; Xu et al, 
2012; Xu et al, 2013a). Analyses of the mitochondrial 
DNA (mtDNA) genome have shown that the tree shrew 
has a closer affinity with Lagomorpha (Arnason et al, 
2002; Xu et al, 2012), while nuclear gene sequences 
provide evidence for a close affinity to primates (Janecka 
et al, 2007; Killian et al, 2001; Lindblad-Toh et al, 2011). 

The tree shrew has long been proposed to be an 
alternative experimental animal to primates in 
biomedical research due to its characteristics such as 
small body size, short reproductive cycle and life span, 
and low-cost maintenance (Peng et al, 1991). This animal 
has been used to create animal models for infection with 
the hepatitis B and C viruses, the development of 
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hepatocellular carcinoma, myopia, psychosocial stress 
(Cao et al, 2003; McBrien & Norton, 1992; van Kampen 
et al, 2002; Yan et al, 1996a; Yan et al, 1996b) and 
learning and memory (Wang et al, 2011), although there 
are many problems awaiting further studies (Xu et al, 
2013b) . Despite an increasing interest in using the tree 
shrew to establish animal models for medical and 
biological research, most of the animals used were 
captured from the wild and/or domesticated for a few 
months. There is no pure inbreeding strain with a clear 
genetic background similar to that of a mouse or rat at 
the present time, which would largely influence the 
stability and repeatability of experiments. 

Tupalldae contains four genera: Tupaia, Anathana, 
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Urogale, and Dendrogale (Helgen, 2005). Currently, a 
total of 15 species is recognized in Tupaia, which is 
broadly distributed across South and Southeast Asia, 
including southern China, India, Philippines, Java, 
Borneo, Sumatra and Bali (Olson et al, 2005; Peng et al, 
1991). Chinese tree shrews (Tupaia belangeri) are 
distributed across southwest China and are divided into 
six subspecies according to geographical distribution and 
morphological characteristics (Wang, 1987). 
Understanding the genetic diversity of the wild tree 
shrew population has critical importance for 
conservation and for a breeding program to establish 
inbred lines. Recently, we evaluated the mtDNA genetic 
diversity of the tree shrew inhabiting urban Kunming and 
observed a relatively high diversity (Chen et al, 2011). 
Genetic markers from the nuclear genome are essential 
for us to further validate this conclusion and to analyze 
other behaviors of this animal, e.g., female-biased 
dispersal and gene flow between different populations. 
Microsatellites are short tandem repeats (STR) in 
the genome and have been broadly used for evaluating 
genetic diversity, population affinity, gene flow and 
population microdifferentiation (Goldstein et al, 1999; 
Schloetterer & Pemberton, 1994; Waits et al, 2000), as 
well as individual identification and parentage tests (Behl 
et al, 2002; Butler, 2005). So far, there have been only 
sporadic reports for characterizing microstatellites in the 
tree shrew (Munshi-South, 2008; Munshi-South & 
Wilkinson, 2006; Srikwan et al, 2002). Amplification of 
the Chinese tree shrew genomic DNA using 
microsatellites developed from other species of 
Scandentia revealed some problems, such as lack of 
polymorphism (He et al, 2009) and had non-specific 
amplification in our samples (data not shown). Therefore, 
it is necessary to develop a set of species-specific 
microsatellite markers based on the available genome of 
the tree shrew. In this study, we successfully retrieved 12 
polymorphic microsatellite markers from Chinese tree 
shrews and characterized them in a population of 117 
individuals. Our results showed that this set of 
microsatellite loci had sufficient data to measure 
population heterozygosity and individual discrimination. 


MATERIALS AND METHODS 


Experimental animals 

A total of 117 Chinese tree shrews inhabiting the 
suburb of Kunming, Yunnan Province, China, were 
collected and raised at the Experimental Animal Core 
Facility of the Kunming Institute of Zoology, Chinese 
Academy of Sciences. All sampling procedures were 
approved by the Institute’s Institutional Animal Care and 
Use Committee. Genomic DNA was extracted from the 
blood, ear and/or muscle tissues of each individual by 
using the genomic DNA extraction kit (Tiangen, Beijing). 
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Microsatellite identification and amplification 

Based on the genomic DNA sequences of the 
northern tree shrew (Tupaia belangeri) provided by the 
Ensembl database (ftp://ftp.ensembl.org/), 20 loci were 
chosen with a core repeat unit of 4-6 nucleotides. Primer 
pairs were designed using primer premier 5. PCR 
amplification condition for each locus contained a pre- 
denaturation cycle at 94 C for 5 min, followed by 30 
cycles of denaturation at 94 ‘C for 30 s, annealing at 
55 C for 30 s, and extension at 72 C for 30 s, followed 
by a final extension cycle at 72 C for 5 min. 
Amplification was performed in a total volume of 20 uL, 
which contained 10xTaq polymerase reaction buffer 
(100 mmol/L Tris-HCl, pH 8.3; 500 mmol/L KCl; 15 
mmol/L MgCl,), 0.2 mmol/L dNTPs, 10 umol/L of each 
forward and reverse primer, 0.5 U Taq DNA polymerase 
(TaKaRa, Dalian, China), and 50 ng genomic DNA. PCR 
products were electrophoresed on 8% polyacrylamide 
gels to test successful and specific amplification. For 
those loci with no amplification or with non-specific 
bands, we optimized the reaction conditions by 
increasing or decreasing the annealing temperature, or 
adjusting the concentration of Mg” to achieve the best 
efficiency and specificity. We validated the authentic 
sequencing of each amplified locus by TA-cloning 
sequencing. In brief, PCR product was ligated into PMD 
19-T vector (TaKaRa, Dalian, China), followed by 
transforming DH5a competent cells (Tiangen, Beijing) 
according to the manufacturer’s instructions. We picked 
up three to five plasmids with the right inserts for 
sequencing. 

After the pretest for primer pairs and the opti- 
mization for amplification, we labeled the forward 
primer of each locus by including 6-FAM (Carboxy- 
fluorescein) at the 5' end of the primer. PCR 
amplification for each locus was performed using the 
above optimized conditions. PCR product was diluted 
with double-distilled water to achieve a concentration of 
1 ng/uL. We mixed 1 uL diluted PCR product with 0.1 
uL Genescan -500 LIZ” (Applied Biosystems) and Hi- 
Di™ Formamide (Applied Biosystems), and added 
ddH2O to make a 10 uL volume of cocktail. The mixture 
was denatured at 95 C for 5 min, then immediately 
chilled on ice. We loaded 1 uL of cocktail on an 
automated sequencer (ABI PRISM 3100, Applied 
Biosystems). Alleles were scored using program 
Genemarker V1.5 (SoftGenetics LLC, State College, PA). 


Statistical analysis 

Allele frequencies, observed heterozygosity (Ho), 
expected heterozygosity (He) and deviation from the 
Hardy-Weinberg equilibrium (HWE) based on the 
likelihood ratio test were calculated by POPgene 1.32 
(Yeh & Boyle, 1997). Polymorphism information content 
(PIC), power of discrimination (PD), power of exclusion 
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(probability of excluding relatives of the true father from 
paternity, PE) and paternity index (PI, which is the ratio 
of the probability that a genetic marker/allele the alleged 
father passed to the child to the probability that a 
randomly selected unrelated man of similar ethnic 
background could pass the allele to the child) were 
calculated by PowerStats V1.2 (Promega, Madison, WI). 
Non-exclusion probabilities (NEP), probability of 
identity (Pid), probability of identity of siblings (Pid-sib) 
and null allele frequency (F (null)) were calculated using 
Cervus 3.0 (Kalinowski et al, 2007). 


RESULTS 


Characterization of the microsatellites 
Among the 20 selected loci, only 14 could be 


successfully amplified in the Chinese tree shrew using 
the primer pairs designed in this study (Table 1). The 
sequence authenticity of each locus was verified by 
sequencing and the microsatellite sequences were 
deposited in GenBank under accession numbers 
JQ173849-JQ173882 (Figure 1). There are some 
mutations in the repeat motif (e.g., from AAG to GAG in 
loci TB9 and TB14), besides the change of the number of 
repeats. The remaining 6 loci either failed to be amplified 
or with apparent non-specific amplification (Table S1). 
To test whether the 14 loci were polymorphic in the 
Chinese tree shrew population, we screened 117 
individuals by using 5'-fluorescently labeled primers. 
Two loci, TB2 and TB10, were non-polymorphic m the 
tested samples. However, a comparison of the sequenced 


Table 1 Information of 14 primer pairs for amplifying microsatellite loci in Chinese tree shrew 
(Tupaia belangeri chinensis) and PCR condition 


Locus Primer sequence (5'-3') 


F:ATCAGAATCTGGTTTCAAAGGT 
TB1 214-224 
R: GCACACCATGATGTAGCTGT 


F:AATATGCTGGTAGGACCAGAAAC 
TB2 185 ° 
R: ATTGGGTTGTTGGGAAAGTC 


F:TCTGGACCTAGCAACTTGCTTTC 
TB3 255-271 
R:TTCTTCACCAGTGTTGAGCATAA 


F:AGACAGAATGCAAGAAATCAC 
TB6 410-431 
R:ATGTGCAATGTAATAGTTCCAG 


F:CTCTGACTTCAAGACCTGGAAC 
TB8 449-473 
R:CTTGAAAGCAAAGGCAGTACAA 


F:CCTTGTGGCTTGAAAAGTATGCT 
TB9 345-389 
R:AGCGAGAATCATCCTCCAGGG 


aT F:ATTGGAGACCACTCTACCCTTTA qe 
R:TGGATTCCCACATAGTCTACCG 


F:GTTTTTCCTCTCTCATGTACTC 
TB12 360-404 
R:GTCCACATCTTTGCATATCTCC 


F:GACCACGGTTCAAGTCCTAG 
TB14 415-479 
R:GATCACAGGGCTGACCTATTA 


F:TGCGGTATCATTGTGGTGGTTCA 
TB15 237-256 
R:GGGCACGTGTGCACATACCTCA 


F:TCCTTAGTATGGACTTGTTTGC 
TB16 206-256 
R:AGGATCACAGGGCAGGGAGACT 


F:ATTACGTCATCCGCTGGTCC 

TB17 290-315 
R:GCAGATTCAGACTTGCCTATCC 
F:TTGACCTAAGAGGTAAATACTG 

TB18 426-498 
R:AAGGTGAGTTTACAGCATAGAC 
F:TGATCTCCAAAGGCTTCTCCAT 


R:GAACCTTCCCGTATGGGTCAGT 


209-229 


Product (bp) 


Annealing temperature ° Mg (mM) Repeat element 

56 °C 1.5 (AAGAAA)n 
‘yee 58°C foradeyeles 15 AAG)M(GAGAAG 
Creo U Clarets j ana 

56 °C l (AAAAC)n...(TTA)n 

55 °C 1 ((A/G)TC(C/T))n 
SERT i woma 

50 °C 1 (ATTT)n 
(68 -60 ‘C)-1°C/cycle for 9 is (T(L/C)(T/C)C )n(PC)n 
cycles, 58 °C for 26 cycles (TTTCCC)n 
cycles 58°C Role 15 AIAG 
ysl 38°C foradeyeles, TOO 
cyeles 58°C foro oyeles 1 AIAG 
mie is 人 > 1 (ATCT)n...(ATATA)n 
myles, 58°C for26cycies 18 TTTC 
(64-60 ‘C)-1°C/cycle for 5 1 (AAAACJn 


cycles, 58 °C for 30 cycles 


“non-polymorphic microsatellite markers in Chinese tree shrews; 
210% DMSO were included in the PCR reaction; 


“PCR amplifications for all loci were performed with 35 cycles. With the exception of TB1, TB6, TB8, and TB10, the other loci were amplified with a 
touch-down protocol. For instance, TB2 was amplified using a decreasing annealing temperature by 1 ‘C every cycle from 64 C to 60 C (total 5 


cycles), then the PCR amplification was performed with an annealing temperature of 58 °C for 30 cycles. 
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Figure 1 GeneScan maps and sequences of 12 polymorphic microsatellites in Chinese tree shrew population 


fragments of these two loci with the reference sequence 
from the Ensembl database showed length 
polymorphisms (Figure S1), suggesting that these two 
loci might present polymorphisms among tree shrews 
from different regions, but not in tree shrews collected 
from Kunming. This result was in accordance with our 
previous studies of mtDNA sequence variation in the 
Chinese tree shrew (Chen et al, 2011; Xu et al, 2012), in 
which we observed some nucleotide differences between 
the Chinese tree shrew from Kunming and the reported 
northern tree shrew outside China, though both were 
grouped as Tupaia belangeri. 


Chinese tree shrew population has a high microsatel- 
lite polymorphism 

Tables 2 and 3 list the estimated parameters for each 
of the 12 polymorphic microsatellites in Chinese tree 
shrews. The number of alleles per locus ranged from 4 to 
20, with an average of 8.75 alleles. The effective alleles 
ranged from 1.19 (TB1 locus) to 10.47 (TB18 locus), and 
had an average of 4.05. The mean values of observed 
heterozygosity and expected heterozygosity were 0.566 
and 0.616, respectively. Loci TBI and TB15 had a 
relatively low level of polymorphism (PIC value<0.25), 
TB3 and TB6 were moderately polymorphic (0.25<PIC 
value<0.5), and the other loci were highly polymorphic 
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(PIC value>0.5) (Table 2). The high PIC value and 
observed heterozygosity indicated a high genetic 
diversity of the analyzed tree shrew population, which 


Table 2 Genetic diversity measures for 12 polymorphic 
microsatellite loci tested in 117 Chinese tree shrews 


Locus Na Ne Ho He PIC 
TB1 4 1.193 0.103 0.107 0.12 
TB3 6 1.786 0.436 0.442 0.40 
TB6 6 2.226 0.537 0.553 0.46 
TB8 7 3.733 0.675 0.735 0.69 
TB9 7 3.256 0.709 0.696 0.65 
TB12 20 7.517 0.607 0.871 0.86 
TB14 17 8.064 0.829 0.880 0.86 
TB15 4 1.213 0.137 0.177 0.17 
TB16 7 3.041 0.675 0.674 0.62 
TB17 5 3.247 0.641 0.695 0.64 
TB18 18 10.47 0.838 0.908 0.90 
TB20 4 2.897 0.607 0.658 0.59 
Average 8.75 4.054 0.566 0.616 0.58 


Na: observed number of alleles; Ne: effective number of alleles; Ho: 
observed heterozygosity; He: expected heterozygosity; PIC: 


polymorphism information content. 
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was consistent with the previous mtDNA analysis (Chen 
et al, 2011). A likelihood ratio test for concordance with 
HWE revealed that only loci TB12, TB17 and TB20 
deviated from HWE (P<0.05), but none was significant 
after applying a Bonferroni correction. 

To further evaluate the potential use of the 12 
microsatellite markers, we calculated some forensic 
efficiency parameters based on the tested samples. The 
linkage disequilibrium was tested by the online SHEsis 
program (Shi & He, 2005), but no significant linkage 
between loci was detected (Table S2). The power of 
discrimination (PD) was from 0.218 (TB1) to 0.971 
(TB18), and the cumulative PD in this study was as high 
as 99.99999998%. Waits et al (2001) recommend a 
reasonable probability of identifying individuals in a 
given region using dominant and codominant markers as 
0.01—0.0001. In this study, the values of combined Pid 
and combined Pid-sib were 4.89x10” and 0.00001, 
respectively. These values are far smaller than 0.01, 
which suggests that these 12 loci had sufficient data for 
individual identification. 

As testing for more markers would increase the cost, 


time and the risk of analytical errors (Waits & Paetkau, 
2005), lt was necessary to choose the minimum number 
of markers that fulfill the aim of individual identification. 
We eliminated the four less polymorphic loci (TB6, TB3, 
TB15 and TB1) but still obtained considerably high 
values for combined PD (CPD=99.99999955%), 
combined Pid (CPid=8.61x10"°), and combined Pid-sib 
(CPid-sib=5.28x 10°). It is thus evident that the remaining 
8 loci are sufficient for individual discrimination. 

For parentage testing, parameters such as average 
non-exclusion probability (NE-IP), average non- 
exclusion probability for one candidate parent given the 
genotype of a known parent of the opposite sex (NE-2P), 
average non-exclusion probability for a candidate parent 
pair (NE-PP), power of exclusion (PE) and paternity 
index (PI) were widely used in the field. The Chinese 
tree shrew population had a NE-1P value of 0.322—0.994, 
a NE-2P value of 0.192—0.945, a NE-PP value of 
0.059—0.897, a PE value of 0.668—0.012, and a PI value 
of 3.05—0.57, respectively (Table 3). Excluding the four 
loci (TB3, TBI5, TBI, and TB6) with few 
polymorphisms, we still obtained sufficient data for 
parentage testing. 


Table 3 Individual identification and parentage testing parameters for 12 polymorphic microsatellite markers in a 
population of Chinese tree shrew (n=117) 


Locus F (null) PI CPI NE-IP CNEIP 


TB18 0.034 3.05 3.05 0.322 0.32 0.192 0.19 0.059 


TB14 0.027 2.93 8.94 0.401 0.13 0.249 4.78x10° 0.095 


TB9 -0.006 1.72 15.37 0.72 9.29x10° 0.544 2.60x10° 0.358 


TB8 0.04 1.58 24.29 0.673 6.25x10° 0.496 1.29x10° 0.309 


TB16 —0.002 4.62x107 7.35x10° 


0.044 


1.54 37.40 0.74 0.57 0.386 


TB17 1.39 51.99 0.727 3.36x107 0.559 4.11x10° 0.38 


TB12 0.19 1.35x107 1.04x10° 


0.046 
0.011 


1.27 66.02 0.404 0.252 0.089 


TB20 1.27 83.85 0.776 1.05x10° 0.621 6.43x107 0.461 


TB6 1.08 90.56 0.845 8.89x10° 0.734 4.72x107 0.599 


TB3 —0.009 7.99x10° 3.60x107 


0.115 


0.89 80.60 0.899 0.763 0.616 


TB15 0.58 46.75 0.984 7.86x10° 0.908 3.27x107 0.83 


TB1 0.048 0.57 26.64 0.994 7.82x10° 0.945 3.09x10” 0.897 


NE-2P CNE2P NE-PP CNEPP 


0.06 


5.54x10° 
1.98x10° 
6.14x10” 
2.37x10” 
9.00x10” 
8.01x10° 
3.69x10° 
2.21x10° 
1.36x10° 
1.13x10° 
1.01x10° 


Pid CPid Pid-sib CPid-sib PD CPD PE CPE 


0.017 0.02 0.302 0.30 0.971 0.971 0.668 0.668 


9.63x10™ 0.966 0.9990140000 
4.22x107 0.849 0.9998511140 
1.74«107 0.886 0.9999830270 


0.028 4.64x107 
0.138 


0.113 


0.319 
0.438 
0.412 


0.654 0.147 


6.43x10° 0.443 0.337 


7.23x10° 0.404 0.605 
7.89x10° 0.837 0.9999972334 


3.48x10° 0.864 0.9999996237 


0.157 1.14x10° 


0.148 


0.454 
0.441 


0.391 0.78 


1.68x107 0.343 0.846 
1.13x107 0.935 0.9999999755 
5.28x10° 0.816 0.9999999955 


2.89x10° 0.71 0.9999999987 


4.67x10° 
8.61x10 7° 


0.028 
0.184 
0.292 


0.032 
0.469 
0.548 


0.299 0.947 


0.299 0.948 


2.52x10!° 0.223 0.968 
1.79x10° 0.653 0.9999999996 


1.49x10° 0.291 0.9999999997 


0.354 8.90x10"! 


0.685 


0.618 
0.833 


0.137 0.979 


6.10x10"! 0.015 0.993 


0.801 4.89x10"' 0.897 1.34x10° 0.218 0.9999999998 0.012 0.995 


F (null): null allele frequency; NE-1P: Average non-exclusion probability for one candidate parent; CNE1P: combined NE-1P; NE-2P: Average non- 
exclusion probability for one candidate parent given the genotype of a known parent of the opposite sex; CNE2P: combined NE-2P; NE-PP: Average 
non-exclusion probability for a candidate parent pair; CNEPP: combined NE-PP; Pid: probability of identity, the probability of that two individual 
share the same genotypes using all loci; CPid: combined Pid; Pid-sib: the probability of identity of siblings; CPid-sib: combined Pid-sib; PD: power 
of discrimination; PE: power of exclusion; CPD: combined PD; PI: paternity index; CPI: combined PI. 


The null allele frequency was also calculated to 
evaluate the efficiency of this set of microsatellite markers 
because of its influence on population analysis (Chapuis 
& Estoup, 2007). Null allele is any microsatellite allele at 
a microsatellite locus that consistently fails to amplify to 
detected levels in the PCR assays, and it has no significant 
effect on parentage analysis when the frequency is less 
than 0.2 (Dakin & Avise, 2004). In our study, all the null 
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allele frequencies were less than 0.2 (Table 3). Therefore, 
these markers can be used for population analysis of tree 
shrews, at least In parentage analysis with little influence 
on the average exclusion probability. 


DISCCUSSION 


Tree shrews have received wide attention in recent 
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years for their potential use in biomedical research and 
drug discovery (Cao et al, 2003; Fuchs & Corbach-Sohle, 
2010; Peng et al, 1991). Despite decades of cultivation 
efforts, no inbred strain has been established for the tree 
Shrew, which hindered a wide use of this animal in 
biomedical research; this might account for the 
inconsistency of previous results concerning animal 
models. Thus establishing some inbred lines of tree 
shrews is crucial, and we have launched an ambitious 
project to fulfill this goal. 

To facilitate our current inbreeding program for the 
Chinese tree shrew and to evaluate the genetic diversity 
of wild tree shrews, we needed a set of microsatellite 
markers that had sufficient data for individual 
discrimination and parentage testing, as well as for the 
quantification of population heterogeneity. In this study, 
by referring to the available tree shrew genome sequence 
(albeit the coverage was only 2x), we were able to obtain 
12 polymorphic loci in the Chinese tree shrew. 
Evaluation of different parameters related to genetic 
diversity, individual identification and parentage testing 
showed that eight of these loci had sufficient data for 
lineage tracing and measurement of diversity. 

Compared with these reported microsatellite 
markers, which were developed for T. glias and T. minor 
(Munshi-South & Wilkinson, 2006; Srikwan et al, 2002), 
our newly designed set based on the genome sequence of 
T. belangeri was specifically optimized for Chinese tree 
shrews. As we had no samples of other species in Tupaia, 
we did not know whether our primer pairs would work 
for these species. As Chinese tree shrews were proposed 
to contain six subspecies (Wang, 1987), it would be 
rewarding to screen the genetic diversity of these 
different subspecies, both for conservation of genetic 
resources and for selecting a proper founder to establish 
the inbred lines. 
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Table S1 Information of the 6 pairs of primers that could not be amplified or had 


non-specific amplification 


Locus 


TB4 


TB5 


TB7 


TB11 


TB13 


Primer sequence (5'-3') 


F:GTCCTAGACCTGGCAGCTTGCT 
R:TTGCTATGTGACCATCCCTGTG 
F:GGGTAGATACCCAGGAGTAGGAT 
R:AGAGATTGCAGTGTGATCATGC 
F:TAAGGACCTGAAAGATGTTTGC 
R:TCCCATGACTGACTGTGGTTTG 
F:TTTGTACTTTATGTAGTTACTC 
R:TTGAAAACACATCAGGCAC 
F:TGGCTTAGACTCCTGCTCGGTG 
R:TTACCACTGCACCATCCTGACG 
F:GTAGCTTGAAAATCATGACAGGC 


R:GAATGCTAGAAACACTCAGTGGC 


Primer 
length 
(bp) 
22 

22 

23 

22 

22 

22 

22 

19 

22 

22 

23 


Tm (°C) 


GC% 


Predicted 
product 
(bp) 

400 

426 

582 


450 


367 


Repeat element 


(AGAT)n 


(CCTTT)n (CCTTC)n 


(TTTC)n (TTTCTC)n 


(TTTC)n(TTCC)n 


(AAAG)n(GAAG)n 


(GAAA)n 
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Table S2 Linkage disequilibrium of 12 polymorphic microsatellite markers in Chinese tree shrew populations 


TBI 
TB3 
TB6 
TB8 
TB9 
TB12 
TB14 
TBI5 
TB16 
TB17 
TB18 
TB20 


TBI 


KK 


0.026 
0.001 
0.012 
0.007 
0.005 
0.006 
0.001 
0.018 
0.022 
0.007 
0.001 


TB3 
0.389 


米 米 
0.017 
0.018 
0.017 
0.006 
0.007 
0.01 
0.02 
0.01 
0.013 
0.008 


TB6 


0.145 
0.217 
KK 
0.007 
0.014 
0.01 
0.01 
0.002 
0.013 
0.009 
0.01 
0.004 


TB8 


0.443 
0.304 
0.17 
米 米 
0.010 
0.007 
0.014 
0.01 
0.012 
0.017 
0.007 
0.017 


TB9 
0.346 
0.308 
0.248 
0.297 
KK 
0.007 
0.008 
0.023 
0.019 
0.008 
0.014 
0.005 


TB12 
0.493 
0.336 
0.393 
0.352 
0.346 
米 米 
0.01 
0.013 
0.009 
0.006 
0.008 
0.008 


TB14 


0.481 
0.289 
0.304 
0.357 
0.3 
0.456 
KK 
0.009 
0.008 
0.009 
0.009 
0.011 


TB15 
0.332 
0.207 
0.134 
0.352 
0.637 
0.519 
0.459 
米 米 

0.021 
0.002 
0.008 
0.006 


TB16 
0.453 
0.254 
0.235 
0.264 
0.302 
0.358 
0.302 
0.484 
KK 

0.003 
0.011 
0.011 


TB17 


0.461 
0.197 
0.176 
0.268 
0.205 
0.306 
0.302 
0.157 
0.122 
KK 

0.012 
0.014 


TB18 


0.572 
0.387 
0.346 
0.349 
0.408 
0.491 
0.47 

0.522 
0.358 
0.368 


KK 


0.012 


TB20 
0.142 
0.17 
0.126 
0.243 
0.119 
0.3 
0.283 
0.236 
0.2 
0.229 
0.349 


米 米 


The upper and lower triangular matrixes present the D values and the r° values, respectively. The values were estimated based on 117 Chinese 


tree shrews. 
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(A) TB2R AATATGCTGGTAGGACCAGAAACAAAGGAAGGAACTAAATAAAGATTTTTTTAAAAGAAG 
TB2C AATATGCTGGTAGGACCAGAAACAAAGGAAGGAACTAAATAAAGATTTTTTTAAAAGAAG 
TB2R AAGAAGAAGAAGAAGAAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAGAAGGAG 
TB2C AAGAAGA--- 一 一 一 一 一 一 一 一 一 一 一 一 AG 
TB2R AAGGAGAAGGAGAAGGAGAAGAAGGAGAAGAAGAGAATGAACATATTAGGTCATCAGGAG 
TB2C AAGGAGAAGGAGAAGGAGAAGAAGGAGAAGAAGAGAATGAACATATTAGGTCATCAGGAG 
TB2R ATTCATGACATATTTTCCTATGCAAAAATGCATCATGACTTTCCCAACAACCCAAT 
TB2C ATTCATGACATATTTTCCTATGCAAAAATGCATCATGACTTTCCCAACAACCCAAT 

(B) TBloR GGACTTTATGCCTCATATCCAACCCTCTTGGGTATTTATTTATTTATTTATTTATTTATT 
TB10C GGACTTTATGCCTCATATCCAACCCTCTTGGGTCTTTATTTATTTATTTATTTATTTATT 
TB10R TATTTATTTATTTATTTATTTATTTATTTTTGGAGAGGTTTTCGCTACAGGGCGCATGCG 
TB10C T 一 -一 一 一 -一 一- 一- TTGGAGAGGTTTTCGCTACAGGGCGCGTGCG 
TB10R CTCGCCAGACAGAGAGACATAGAGACAGAAAGCGAGAAATCACCCCCCAGAGGGGATGGT 
TB10C CTCGCCAGACAGAGAGACATAGAGACAGAAAGCGAGAAATCACCCCCCAGAGGGGATGGT 
TBLOR CACGCCGCAATTCCTICTCCTCCATICTCCTCTCGGGATCTTGGGGACGCACGCACGCGC 
TBLOC CACGCCGCAATTCCTTCTCCTCCATTCTCCTCTCGGGATCTTGGGGACGCACGCACGCGC 


Figure S1. Alignment of sequences of loci TB2 (A) and TB10 (B) in Chinese tree 
shrew with the tree shrew reference from the Ensembl database. There are differences 
for the number of repeat motif in Chinese tree shrews (locus TB2: TB2C; TB10: 
TB10C) and the tree shrew from the Ensembl database. 
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Abstract: Establishing non-human primate models of human diseases is an efficient way to narrow the large gap between basic 
studies and translational medicine. Multifold advantages such as simplicity of breeding, low cost of feeding and facility of operating 
make the tree shrew an ideal non-human primate model proxy. Additional features like vulnerability to stress and spontaneous 
diabetic characteristics also indicate that the tree shrew could be a potential new animal model of human diseases. However, basal 
physiological indexes of tree shrew, especially those related to human disease, have not been systematically reported. Accordingly, 
we established important basal physiological indexes of domesticated tree shrews including several factors: (1) body weight, (2) core 
body temperature and rhythm, (3) diet metabolism, (4) locomotor rhythm, (5) electroencephalogram, (6) glycometabolism and (7) 
serum and urinary hormone level and urinary cortisol rhythm. We compared the physiological parameters of domesticated tree shrew 
with that of rats and macaques. Results showed that (a) the core body temperature of the tree shrew was 39.59+0.05 °C, which was 
higher than that of rats and macaques; (b) Compared with wild tree shrews, with two activity peaks, domesticated tree shrews had 
only one activity peak from 17:30 to 19:30; (c) Compared with rats, tree shrews had poor carbohydrate metabolism ability; and (d) 
Urinary cortisol rhythm indicated there were two peaks at 8:00 and 17:00 in domesticated tree shrews, which matched activity peaks 
in wild tree shrews. These results provided basal physiological indexes for domesticated tree shrews and laid an important foundation 
for diabetes and stress-related disease models established on tree shrews. 


Keywords: Domesticated tree shrew; Basal physiological indexes; Rhythm; Glycometabolism; Stress-related hormone 


With fast proliferation and low cost, it is also a potential 
non-human primate model proxy. Furthermore, many of 


Translational medicine always faces the dilemma of 
gaps between experimental rodents and humans. Non- 
human primates are considered the best animal models 
for human diseases. However, macaque models are 


largely limited by expense, long experimental cycle, Received: 11 December 2012; Accepted: 06 February 2013 


expert shortages, and strict animal rights. Thus, it is 
necessary to find an economical and easily operated non- 
human primate animal model. 

The tree shrew (Tupaia belangeri chinensis) 1s 
considered a close primate sister, with high genetic 


similarity to humans (Elliot O, 1971; Janecka et al, 2007). 
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its features make the tree shrew a suitable model for 
human disease, especially in relation to stress-related 
diseases (Cao et al, 2003; Carter, 2007). For example, 
chronic social defeated depression is well established for 
tree shrew due to its vulnerability to stress (Fuchs, 2005; 
Wang et al, 2011). In addition, the existence of spontaneous 
diabetes in tree shrews makes it possible to establish 
diabetic models for this species (Liang & Wu, 2011). 

Compared with well-investigated experimental 
animals, the physiological parameters of T. b. chinensis 
have been rarely studied. Here we tested a set of basal 
physiological indexes for the tree shrew, including daily 
food and water intake, systolic and diastolic blood 
pressure, fasting blood glucose, and basal levels of 
testosterone, estradiol and cortisol. Understanding these 
parameters could be useful for a variety of research 
studies using the tree shrew as the model system. 


MATERIALS AND METHODS 


Animals 

We used adult tree shrews (n=92, 130-150 g) 
obtained from a breeding colony at the Animal House 
Center of the Kunming Institute of Zoology, and 
Sprague-Dawley (SD) rats (n=8, 250-300 g) obtained 
from the Kunming Medical University. Animals were 
individually housed under a 12 h/12 h dark/light cycle 
(light, 8:00-20:00; dark, 20:00—8:00) in thermoregulated 
rooms (T: 25-27 °C, RH: 55%-—70%) with free access to 
water and food. Animal care and experimental protocols 
were approved by the Animal Ethics Committee of the 
Kunming Institute of Zoology, the Chinese Academy of 
Sciences. 


Body-mass and cardiovascular index 
Body weight, 24-hour diet, locomotion and temperature 
rhythm 

Eight tree shrews were weighed in the morning. 
Their food and water intake were measured In metabolism 
cages. In consideration of food moisture evaporation 
factors, we also placed a plate of food outside the cage as 
a control. After six days of 24-hour diet measurement, 
the tree shrews were returned to their home cage (size: 
100 cm x 68 cm x 86 cm, [wxdxh]) consisting of a glass 
side. Subjects were videotaped from 8:00 to 20:00 and 
their movement distance in the vertical plane was 
recorded by a Noldus EthoVision XT Version 8.0 video 
tracking system (Wageningen, the Netherlands). Core 
body temperature was measured using an anal 
thermometer probe, inserted 2 cm from the anal sphincter 
of the tree shrew subjects. 


Telemetry monitoring of blood pressure and heart 
rate 


Blood pressure and heart rate of tree shrews were 
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tested by telemetry monitoring, as described previously 
(Coolen et al, 2012). In brief, adult tree shrews (n=4, 3 
male and 1 female) were anaesthetized with 
pentobarbital sodium (100 mg/kg, 1.p.). Transmitters with 
radiotelemetric catheters were implanted in enterocoelia 
to measure blood pressure and heart rate (TAC50-PXT; 
Data Sciences International, St Paul, MN, U.S.A.). One 
radiotelemetric catheter was inserted into the abdominal 
aorta below the renal arteries and fixed into place with 
tissue adhesive. The implant was then sutured to the 
abdominal wall, and the abdominal wall and skin 
incision were sutured with silk suturing. When the 
animals were awake, blood pressure and heart rate were 
detected by a PhysiolTel-Receiver (MOPEL RPC-1; Data 
Sciences International, St Paul, MN, U.S.A.) installed 
under each individual sleeping-box. Data signals were 
analyzed using Dataquest A.R.T. t Gold Acquisition, 
Version 2.20. 


Stress-related hormone and sexual hormone 
Urine cortisol samples 

Urinary cortisol rhythm was tested in tree shrews 
(n=30). Urine samples were collected each hour from 
08:00 to 20:00 and kept at —20 °C. Samples were 
measured with an Iodine ['”°I] cortisol radioimmunoassay 
kit (Beijing North Institute of Biological Technology, 
China) and y radio-immunity counter (GC-2010, Zonkia, 
Anhui, China), as per manufacturer’s instructions. 


Enzyme-Linked Immuno Sorbent Assay (ELISA) 

Whole blood (about 0.2 mL) was obtained from the 
femoral vein of tree shrews (n=28) using a 2 mL syringe, 
and centrifuged at 4 C for 8 min at 5 000 r/min. The 
supernatant was taken and stored at —20 °C. Sex 
hormones such as testosterone (T) and estradiol (E2), and 
serum level of stress-related hormones such as cortisol, 
norepinephrine (NE), adrenocorticotrophic hormone 
(ACTH) and corticotrophin releasing hormone (CRH) 
were assayed using human ELISA kit (Shanghai Li-Chen, 
China), as per manufacturer’s protocols. 


Glycometabolism index 
Blood glucose and oral glucose tolerance (OGTT) 
Blood glucose was tested from tail blood using a 
ONE TOUCH BASIC meter (Lifescan, Inc 2001 
Milpitas, CA 95035, USA). Eight tree shrews and eight 
SD rats were fasted for 12 h prior to the experiment. 
Thereafter, fasting blood glucose (FBG) and 1-hour 
postprandial blood glucose (PBG) after feeding were 
tested. In OGTT, animals were given glucose (2 g/kg, 4 
g/kg) or saline (1 mL/kg) by intragastric administration 
(i.g.) after 0 point test. Blood glucose was detected at 30 
min, 60 min, 90 min, 120 min, and 360 min after 
intragastric treatment. 
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Sucrose consumption 

Tree shrews were individually housed (n=8) and 
were given two bottles containing 200 mL of water and 
200 mL of sucrose solution. Each concentration (0.5%, 
1%, 2%, 5%, 10% and 20%) was measured once a day, 
for 2 days. The volumes of consumed water and sucrose 
solution were recorded. 


Electroencephalograms (EEG) 

Three male and three female tree shrews were 
anaesthetized with pentobarbital sodium (100 mg/kg, 
1.p.), ventilated with 95% O, and 5% CO, and with their 
core temperature maintained at 3740.5 °C. Animals were 
stereotaxically (RWD Life Science Co., China) 
implanted with recording electrodes in the hippocampus 
(15.6 mm posterior to bregma, 6 mm lateral to midline, 
9.5 mm below dura) and cerebral cortex (15.6 mm 
posterior to bregma, 6 mm lateral to midline, 2 mm 
below dura) on the right hemisphere of the tree shrew 
brain to measure EEG under anesthetic. Recording 
electrodes were made by gluing together a pair of twisted 
Teflon-coated 90% platinum and 10% iridium wires (50 
um inner diameter, 75 um outer diameter; World 
Precision Instruments, Sarasota, FL). The EEG signals 
were amplified and filtered as follows: high-pass filter at 
0.1 Hz and low-pass filter at 2 kHz. Sample frequency 
was 200 Hz. The EEG power spectra were computed by 
a Fast Fourier Transform routine for 4-s epochs (0.25 Hz 


weight, average food and water intake per day, basal 
levels of T, E2, and cortisol were tested. We also 
measured the systolic and diastolic blood pressure, heart 
rate and core body temperature by telemetry monitoring 
in freely moving tree shrews (Table 1). Compared with 
published data of other well established laboratory 
animals, such as rats and macaques, the physiological 
parameters in our study showed that tree shrews had the 
highest heart rate and core body temperature, which was 
likely related to its high activity level. Systolic pressure 
of tree shrews was much closer to macaques. Also 
similar to macaques, the main stress-related hormone of 
tree shrew was cortisol, compared to corticosterone in 
rats (Table 2). 


Table 1 Basal parameters of domestic tree shrews 


Indexes Tree shrew 
Body-mass index Body weight (g) 147.90+1.23 
Food intake (g/d) 116.40+8.09 
Water volume (mL/d) 4.33+41.38 
Core body temperature (°C) 39.59+0.05 
High activity time 17:30 — 19:30 
Cardiovascular Systolic pressure (kPa) 20.92+0.04 
Diastolic pressure (kPa) 11.8142.88 


Heart rate (beats/min) 357.0043 1.00 


resolution). Sexual hormone Testosterone (T, ng/mL) 0.1340.01 
Estradiol (E2, ng/mL) 42.29+3 .36 
Data analysis Gl tabolism Fasting blood gl (FBG, W/L) 6.20+0.84 
ee ` cometabolism Fastin OO UCOSe mmo š Fi 
Data were statistically analyzed using SPSS 19.0. 和 
Blood glucose data were analyzed by a Paired-Sample t- Posprandialiplood gu se (EDC mn sto do 
test and sugar preference concentration test data were Sugar sensitive concentration 5% 
analyzed by One-Way ANOVA test, followed by post- R dolom 200.8941424 
hoc comparisons (Fisher’s LSD test). Significance level a a 
y 
was set at P<0.05. Data was expressed as mean+SEM. Morepinepunne (NE aon) HRN 
Adrenocorticotrophic hormone 
(ACTH, pg/mL) §7.3543.17 
RESULTS Corticotrophin releasing hormone 
64.55+2.13 
B ] boli ind (CRH, pg/mL) 
asal metabolism indexes 
: i Urinary cortisol (ng/mL) 0.63+0.06 
Basal metabolism index measurements, body 
Table 2 Basal parameters of tree shrew compared with rats and macaques 
Indexes Rat Tree shrew Macaque 
Body-mass index Core body temperature (°C) 38.80+40.21 39.59+0.05 38.70+40.26 
Cardiovascular Systolic pressure (kPa) 17.29+1.06 20.92+0.04 24.1043 .67 
Diastolic pressure (kPa) 11.70+2.13 11.8142.88 18.7943 .50 
Heart rate (beats/min) 289.86+18.36 357.00+31.00 157.92+12.59 
Glycometabolism FBG (mmol/L) 7.80+2.67 6.20+0.84 4.65+40.53 
Stress-related hormone Cortisol (ng/mL) 200.89+14.24 107.80+17.60 
Corticosterone (ng/mL) 171.18+439.98 = = 
ACTH (pg/mL) 22.88+7.41 57.3543.17 77.88+21.34 


Data on rats were from Fan et al (2012), Li et al (2011), Liu et al (2006), Xu et al (2006) and Yu et al (2006), and data on macaques were from Xu et al (2007) 


and Zhang et al (1998). 
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Electroencephalograms (EEG) 

The EEG power spectra of anaesthetized tree shrews 
were recorded by an electrophysiological recording 
system (Figure 1A), which showed there were no gender 
differences in the EEG of domesticated tree shrews. 
Delta waves (6, 1—4 Hz) were the major waveform in the 
hippocampus and cerebral cortex (Figure 1B, C). 


Circadian rhythm 

For circadian rhythm measurement, tree shrews 
were videotaped from 8:00 to 20:00 and their movement 
distance in the vertical plane was recorded. Results 


IFA 


showed that tree shrews adapted their behavior in an 
uncertain environment, with animals fed in the 
laboratory showing an activity peak period only from 
17:30 to 19:30 (Figure 2A). Urine cortisol rhythm was 
determined by urine sampling at 1-hour intervals from 
08:00 to 20:00. Results showed there were two activity 
peaks at 8:00 and 17:00, respectively, which matched 
activity peaks in the wild (Figure 2B). Core body 
temperature was recorded by anal temperature method 
similar to that often used in rats. During the daytime, 
minimum body temperature was about 38.7 C and 
maximum was about 40 °C (Figure 2C). 


Hippocampus. 


Cerébral cortex 
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Figure 1 Electroencephalogram of the hippocampus and cerebral cortex in tree shrews 
EEG power spectra showed no gender differences in tree shrews. The 5 waves were the major waveform of the hippocampus (B) and cerebral cortex (C) in 


anaesthetized tree shrews. 
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Figure 2 Circadian rhythm of locomotion, urinary cortisol and core body temperature 
A: Locomotion activity peak period was from 17:30 to 19:30; B: Urinary cortisol showed two peaks at 8:00 and 17:00 respectively; C: Core body temperature 


showed a minimum body temperature of 38.7°C and maximum of 40°C. 
Glycometabolism 


Blood glucose was tested from tail blood of tree 
shrews and l-hour PBG after feeding was double the 
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FBG (Figure 3A, #=—3.442, P=0.014). The OGTT results 
showed that given 2 g/kg glucose to fasted animals, the 
blood glucose of tree shrews increased rapidly for the 
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Figure 3 Detection of blood glucose and sugar concentration for tree shrews 
A: PBG was significantly higher than FBG; B: Tree shrews and SD rats were given glucose or saline (i.g.) after 0 point test in OGTT. Blood glucose was tested 
at 30, 60, 90, 120 and 360 min after intragastric treatment; C: 5% was the best sugar concentration among the 0.5, 1, 2,5, 10 and 20%. *:P<0.05, **:P<0.01. 


first 30 min. After 6 hours metabolizing, the blood 
glucose level was still high (Figure 3B). The sugar 
preference concentration test indicated that tree shrews 
preferred the 5% sugar concentration (Figure 3C, 
Fis 885.964, P=0.000). 


DISCUSSION 


Due to the large differences between rodents and 
humans, it is extremely difficult to translate basic 
experimental findings into clinical treatment. Thus, we 
need new animal models that can reveal mechanisms 
more closely related to those of human disease. Tree 
shrews are closely related to primates, and are 
characterized by fast proliferation, short experiment 
cycle, and low cost. The brain functions of tree shrews 
are also very similar to those of humans, and established 
standardized tree shrew strains could be used to study 
human disease (Wang et al, 2012). Considering the 
limited studies previously conducted on T b. chinensis, 
we determined certain basal physiological indexes for 
domesticated tree shrews, including body weight, 24- 
hour diet, locomotion and temperature rhythm, systolic 
and diastolic blood pressure, heart rate, blood glucose, 
OGTT, sucrose consumption (0.5%, 1%, 2%, 5%, 10% 
and 20%), basal levels of sex hormones such as 
testosterone and estradiol and stress-related hormones 
such as cortisol, norepinephrine, ACTH and CRH. The 
EEG of the hippocampus and cerebral cortex were also 
measured. 

Compared with rats and macaques, physiological 
parameters showed that tree shrews had the highest heart 
rate and core body temperature. Results demonstrated 
that tree shrews had a high activity level, which makes it 
an ideal animal model for metabolism studies. Similar to 
macaques, the main stress-related hormone of tree 


shrews was cortisol, compared with corticosterone in rats. 


Circadian rhythm measurement showed that both tree 
shrews and macaques were diurnal, while the rat was 
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nocturnal. These characteristics suggest that the tree 
shrew is closely related to non-human primates and is a 
good animal model for stress-related disease. 

For circadian rhythm measurement, domesticated 
and wild tree shrews had similar cortisol rhythm peaks 
but different activity peak times. That is, the activity 
peak period for tree shrews fed in the laboratory showed 
only one peak from 17:30 to 19:30, but did not show the 
second peak between 8:00 to 10:00 as seen in wild tree 
shrews (Fuchs & Corbach-Sohle, 2010). However, 
cortisol rhythm in the domestic tree shrew showed two 
peaks at 8:00 and 17:00, consistent with activity peak 
times in the wild. The main reason for this difference 
was because the breeders cleaned the animal house from 
8:00 to 10:00, and thus the domesticated tree shrew had 
adapted to the laboratory environment. 

For sucrose consumption, the 5% concentration was 
favored by the tree shrews, compared with the 1% 
sucrose concentration preferred by rodents (Tye et al, 
2013). In the oral glucose tolerance test, tree shrews 
showed more sensitivity to glucose than rats but 
demonstrated poor carbohydrate metabolism ability. 
These might be potential causes for the exhibition of 
spontaneous diabetes in tree shrews (Liang & Wu, 2011). 
Tree shrews would be a good animal model for 
investigating etiology and pathophysiology of diabetes in 
future studies. 

Our results broaden understanding of the 
physiological characteristics in tree shrews and also 
provide detailed reference for future basic scientific 
research and preclinical study based on tree shrews. 
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中 国 科 学 院 昆 明 动 物 研 究 所 遗传 资源 与 进化 国家 重点 实验 军 
2013 年 “进化 生物 学 ” 嗜 期 班 遂 知 


中 国 科 学 院 昆明 动物 研究 所 遗传 资源 与 进化 国家 重点 实验 室 将 于 2013 年 8 月 9 日 一 18 日 举办 “进化 生 
物 学 ” 香 期 班 。 

晋 期 班期 间 将 邀请 本 领域 著名 专家 学 者 进行 进化 生物 学 专题 讲座 并 邀请 本 所 专家 介绍 本 领域 研究 前 
治 进 展 ， 组 织 学 员 开 展 经 典 论 文 的 报告 和 辩论 ， 促 进 学 员 间 的 学 习 交 流 ， 引 寻 和 培养 学 生 思 考生 物 学 问题 
的 正确 方法 ， 还 将 组 织 学 生 参 观 址 传 资 源 与 进化 国家 重点 实验 室 、 昆 明 动 物 博 物 馆 等 。 

进化 生物 学 一 直 是 生命 科学 领域 的 核心 问题 之 一 。 近 年 来 ， 随 着 基因 组 学 的 迅 独 发 展 ， 进 化 生物 学 得 
到 了 突破 式 发 展 。 中 国 科 学 院 昆 明 动 物 研 究 所 遗传 资源 与 进化 国家 和 章 点 实验 室 在 遗传 资源 多 样 性 的 演化 与 
保护 、 基 因 与 基因 组 的 进化 、 遗 传 发 育 与 进化 等 方面 的 部 分 研究 处 于 国际 领先 地 位 。 目 2002 年 起 ， 实 验 
室 已 举办 了 9 期 “昆明 进化 生物 学 嗜 期 学 习 研 讨 班 >， 得 到 了 广大 学 生 和 相关 研究 人 员 的 普 过 赞誉 。 全 今 共 
WU ia 745 人 ， 其 中 四 分 之 三 来 目 985 院 校 。 目 前 ， 已 经 有 不 少 嗜 期 班 学 员 在 国内 、 外 高 校 中 成 为 优秀 
的 青年 科学 家 ， 包 括 两 名 “青年 和 后 人 ”和 多 名 正教 授 。 

本 次 吐 期 班 由 中 国 科 学 院 大 学 和 遗传 资源 与 进化 国家 重点 实验 室 共同 资助 。 我 们 将 在 全 国 主要 高 校 戏 
选 45 名 大 二 以 上 的 本 科 生 参加 ， 对 正式 入 选 的 同学 ， 我 们 将 报销 往返 硬座 和 车票， 并 提供 在 昆明 学 习 期 间 
HIRIE BY o 

申请 的 本 科 生 同学 必须 热爱 生命 科学 ， 对 有 关 学 科 有 浓厚 兴趣 ， 并 提供 以 下 材料 : 
1. 个 人 简历 一 份 ， 介 绍 目 己 的 高 考 成 绩 、 专 业 学 习 情 况 、 诛 外 阅读 和 社会 活动 、 获 奖 情况 、 通 信和 方 
2. 提交 一 封 个 人 申请 说 明 (personal statement) 和 一 篇 关于 生物 和 学、 基因 组 学 或 进化 生物 学 的 综述 
或 讨论 性 论文 ， 字数 不 限 。 
3. 本 次 赫 期 班 将 采取 科研 实践 与 授 读 相 结合 的 形式 ， 参 加 学 生 需 根据 个 人 兴趣 与 重点 实验 室 各 诛 题 
组 研究 方向 ， 选 择 1 一 2 个 读 题 组 作为 亚 期 班期 间 主 要 的 科研 实践 实验 室 ， 并 参与 实验 室 的 科研 
活动 .请 有 意 参 加 的 同学 提前 与 各 读 题 组 负责 人 联系 。 实 验 军 相关 信息 请 见 国 家 重点 实验 室 网 站 : 


http://www.kiz.cas.cn/gre/ 








































































































请 将 上 述 材料 以 电子 邮件 方式 于 6 30 日 前 寄 到 会 议 联系 人 信箱 。 我 们 将 认真 评审 ， 于 7 月 10 日 前 
通知 入 选 的 同学 。 


欢迎 有 意 报考 本 实验 室 研究 生 或 有 免 推 资格 的 同学 积极 参加 ! 请 对 此 感 兴趣 的 同学 定期 关注 本 所 香 期 
班 网 站 的 更 狐 信 息 。 哮 期 班 网 站 :http:/www.kiz.cas.cn/gre/sSummercourse/ 


联系 人: 唐 8 
电子 邮件 : tangj@mailkiz.ac.cn 
电 话 : 0871-65199125 


地 址 : 云南 省 昆明 市 教 场 东 路 32 号 
HB 4: 650223 





